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In the free three-dimensional space, we consider a pair of identical t fermions of some species 
or in some internal state, and a pair of identical fermions of another species or in another state. 

There is a resonant s-wave interaction (that is of zero range and inhnite scattering length) between 
fermions in different pairs, and no interaction within the same pair. We study whether this 2 -|- 2 
fermionic system can exhibit (as the 3-1-1 fermionic system) a four-body Efimov effect in the absence 
of three-body Ehmov effect, that is the mass ratio a between d and fermions and its inverse are 
both smaller than 13.6069. ... For this purpose, we investigate scale invariant zero-energy solutions 
of the four-body Schrodinger equation, that is positively homogeneous functions of the coordinates 
of degree s — 7/2, where s is a generalized Efimov exponent that becomes purely imaginary in the 
presence of a four-body Efimov effect. Using rotational invariance in momentum space, it is found 
that the allowed values of s are such that M(s) has a zero eigenvalue; here the operator M{s), that 
depends on the total angular momentum i, acts on functions of two real variables (the cosine of 
the angle between two wave vectors and the logarithm of the ratio of their moduli), and we write 
it explicitly in terms of an integral matrix kernel. We have performed a spectral analysis of M{s), 
analytical and for an arbitrary imaginary s for the continuous spectrum, numerical and limited to 
s = 0 and £ < 12 for the discrete spectrum. We conclude that no eigenvalue of M(0) crosses zero 
over the mass ratio interval a G [1,13.6069 . ..], even if, in the parity sector (—1)^, the continuous 
spectrum of M(s) has everywhere a zero lower border. As a consequence, there is no possibility of 
a four-body Ehmov effect for the 2-1-2 fermions. 

We also enunciated a conjecture for the fourth virial coefhcient of the unitary spin-1/2 Fermi gas, 
inspired from the known analytical form of the third cluster coefHcient and involving the integral 
over the imaginary s-axis of s times the logarithmic derivative of the determinant of M{s) summed 
over all angular momenta. The conjectured value is in contradiction with the experimental results. 

PACS numbers: 67.85.-d, 21.45.-v, 34.50.-s 


I. INTRODUCTION AND MOTIVATION 

In three-dimensional cold atomic gases, thanks to mag¬ 
netic Feshbach resonances, it is now possible to induce 
resonant s-wave interactions between the particles [l|. 
This means that the s-wave scattering length a is in ab¬ 
solute value much larger than the range (or the effective 
range) of the interaction. Essentially, one can assume 
that 1/a = 0, and since the de Broglie atomic wavelength 
is also much larger than the range of the interaction, one 
can replace the interactions by scaling invariant Wigner- 
Bethe-Peierls two-body contact conditions on the wave- 
function [^: one realizes the long sought unitary limit. 

Perhaps the most striking phenomenon that can take 
place in that regime is the Efimov effect, predicted for 
three particles with appropriate statistics and mass ra¬ 
tios 3 . It corresponds to the occurrence of an infinite 
number of bound states, with an asymptotically geomet¬ 
ric spectrum close to the zero-energy accumulation point. 
The geometric part of the spectrum is characterized by 
a ratio, predicted by Efimov’s zero-range theory, and a 
global energy scale that depends on the microscopic de¬ 
tails of the interaction. The mere existence of such an 
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energy scale forces us to supplement the two-body con¬ 
tact conditions by three-body ones, that involve a length 
scale, the so-called three-body parameter, and that break 
the scale invariance at the three-body level. It is at this 
cost that the zero-range model becomes well defined and 
leads to a self-adjoint Hamiltonian. The Efimov effect is 
now observed experimentally with cold atoms Q, which 
gives access to the value of the three-body parameter @ . 

A natural question is to know whether or not a four- 
body Efimov effect is possible leading to an in¬ 

finite, asymptotically geometric, spectrum of tetramers, 
with an energy ratio predicted by a zero-range theory 
and a global energy scale fixed by a four-body parame¬ 
ter appearing in four-body contact conditions. It is now 
understood that a prerequisite to the four-body Efimov 
effect is the absence of three-body Efimov effect: it is in¬ 
deed expected that the introduction of three-body con¬ 
tact conditions (in terms of the three-body parameter) 
imposed by the three-body Efimov effect is sufficient to 
also render the four-body problem well defined, that is 
without the need for a four-body parameter; as predicted 
in reference Q , no geometric sequence of tetramer states 
can then be found but, as shown numerically for four 
bosons 0, sequences of four-body complex energy reso¬ 
nances are expected in general, with the sa me g eometric 
ratio as the trimer Efimov spectrum (see [T^ for early 
studies not accessing the imaginary part of the energy). 
This prerequisite rules out systems with more than one 
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boson 3 as possible candidates for a four-body Efimov 
effect, and suggests to use fermions to counterbalance the 
Efimov effect by the Pauli exclusion principle, at least in 
three dimensions (what happens in lower dimensions or 
with resonant interactions in other channels than the s- 
wave is discussed in 

Consider then the so-called p + q fermionic problem: p 
identical fermions of the same species or spin state reso¬ 
nantly interact in free space with q identical fermions of 
another species or spin state. It is assumed that there is 
no interaction between the identical fermions, since they 
cannot scatter in the s wave. It is convenient to adopt 
a pseudo-spin notation, with f for the first species and 
I for the second. The two species have in general dif¬ 
ferent masses and m^, and the crucial idea is to use 
their mass ratio as an adjustable parameter to search for 
the four-body Ehmov effect without triggering the three- 
body one. 

The 3-1-1 or ttti resonant fermionic problem was inves¬ 
tigated in reference [l^. A four-body Efimov effect was 
predicted for a mass ratio 13.384 < m-fjmi < 13.6069 .... 
Beyond 13.6069 ... the three-body Efimov effect sets in as 
shown in iQIil, which blocks the four-body Efimov 
effect as discussed above: apart from a finite number of 
tetramer states, one expects an infinite number of four- 
body resonances with the same geometric ratio as for the 
2-1-1 problem. 

The main motivation of the present work is to deter¬ 
mine the presence or the absence of a four-body Efimov 
effect in the 2 -|- 2 or ttii fermionic problem. To our 
knowledge, no general and rigorous answer was given to 
this problem. One may think attacking it with the Born- 
Oppenheimer approximation. We indeed expect (as for 
the three-body case) that the only possibility for a four- 
body Efimov effect is to have a large mass imbalance 
between the two species, for example the f fermions are 
much heavier than the j, ones. It is found that, in pres¬ 
ence of two t fermions at hxed positions, there is a single 
bound state for the j, particle, which creates an effec¬ 
tive oc attraction between the f fermions. 

For a large enough m^/m^ mass ratio, this indeed beats 
the centrifugal barrier oc f{}/{rn^r^) between the f par¬ 
ticles (they are fermions and approach each other with 
a non-zero angular momentum), which qualitatively ex¬ 
plains the occurrence of a three-body Efimov effect in 
the 2-1-1 problem, as pointed out in 1973 by Efimov Q. 
However, as there is a single bonding orbital, one cannot 
put a second ), fermion in that orbital, but one can at 
best put one in the ground, zero-energy scattering state, 
which has two consequences: (i) the Born-Oppenheimer 
attractive potential between the f particles is not low¬ 
ered by the second fermion, so no four-body Efimov 
effect is predicted at a mass ratio strictly below the three- 
body Efimov effect threshold, and (ii) as emphasized in 
the second fermion, being in a zero-energy eigen¬ 
state, does not have a fast motion as compared to the 
one of the heavy particles, which sheds doubts on the 
validity of the Born-Oppenheimer approximation. Alter¬ 


natively, one may expect that this 2 -|- 2 problem was 
already solved numerically in the literature; however no 
convincingly dense coverage of the mass ratio interval 
between 1 and 13.6069... seems to be available in the 
numerics considering the narrowness of the above 
mentioned mass interval. To obtain a firm answer to the 
question, we generalize the method of reference Bin, 
deriving from the zero-range model momentum space in¬ 
tegral equations for the 2-1-2 fermionic problem at zero 
energy (see also the most general formulation of reference 
Ml and using rotational symmetry and scale invariance 
to reduce them to a numerically tractable form. 

Another motivation is to pave the way for the calcu¬ 
lation of the fourth virial coefficient of a two-component 
unitary Fermi gas: this would make an interesting bridge 
between few-body and many-body physics. For a unit 
mass ratio m^/m^ = 1, the value of this virial coeffi¬ 
cient was already obtained experimentally from a mea¬ 
surement of the equation of state of a gas of ultracold 
atoms BHH. On the theory side, there exist two main 
techniques. First, the diagrammatic technique, used ex¬ 
actly (ffil diagrams are kept) for the third virial coeffi¬ 
cient [22Ll23j| . and approximately (only some diagrams are 
kept, those relevant in the perturbative regime of a large 
effective range or a small scattering length) for the fourth 
virial coefficient [2^ leading to a value different from but 
reasonably close to the experimental value. Second, the 
harmonic regulator techni que 12511 . used with success for 
the third virial coefficient [261 - 1^ , that requires to deter¬ 
mine the spectrum of up to four particles in an isotropic 
harmonic trap. A first, brute force numerical solution of 
this trapped four-body problem [3^ was not able to re¬ 
cover even the sign of the experimental value. In a more 
analytical way, this spectrum can be deduced from the 
solutions of the zero-energy free space problem [3l|, , 

due to the SO(^l) dynamical symmetry of the unitary 
Fermi gas [33 - [^ . so that the four-body integral equa¬ 
tions written here may also be useful for the solution of 
the virial problem. 

Our article is organized as follows. In section [II] we 
derive the zero-energy momentum-space integral equa¬ 
tions in general form. In section IIIII we successively use 
the rotational invariance, the scale invariance and the 
parity invariance to put the integral equations in a max¬ 
imally reduced form. This reduced form, written in sec¬ 
tion [IV^ exactly expresses the fact that some operator 
M, depending on the angular momentum ^ and the scal¬ 
ing exponent s, has a zero eigenvalue, which motivates its 
spectral analysis; it allows to show that two components 
of the continuous spectrum of M can be expressed exactly 
in terms of the Ehmov transcendental functions appear¬ 
ing in the fti and fii three-body problems (see section 
IIV Bl) and that there is a third, unexpected continuum 
due to a term with no equivalent in the 3-1-1 problem 
(see section IIV Cl) . The question of the existence of the 
four-body Ehmov effect in the 2 -(- 2 fermionic problem 
is the subject of section El whereas the secondary mo¬ 
tivation of this work, i.e. the fourth virial coefficient of 
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the spin 1/2 unitary Fermi gas, is relegated to the Ap¬ 
pendix |51 where its expression in terms of the operator 
M is conjectured from a transposition of the known ana¬ 
lytic expression of the third virial coefficient [ 2 ^, [ 2 ^ , and 
the conjectured value is compared to the experimental 
[13, [m and theoretical [13, values. We conclude in 
section ED 


II. DERIVATION OF THE GENERAL 
FOUR-BODY INTEGRAL EQUATIONS 

Particles 1 and 2, of positions ri and r 2 , belong to 
species t- Particles 3 and 4, of positions ra and r 4 , 
belong to species j,. The four-body wavefunction ip is 
subjected to the usual Wigner-Bethe-Peierls contact con¬ 
ditions, for a zero-range interaction of s-wave scatter¬ 
ing length a between opposite-spin particles. For all 
i G {1,2} and all j € {3,4}, when the distance be¬ 
tween particles i and j tends to zero, at fixed position 
Rij = (rrifri + m_i^rj)/(m-f + mi) of their center of mass 
(different from the positions of the remaining two parti¬ 
cles), one imposes 


r-ij-s-o \rij a) zirm 
^ "^ij)k^i,j) 4" ^i.^ij) ( 1 ) 

where the form of the regular part Aij supposes that 
the center of mass of the four particles is at rest, and 
where = mimi/(mi+mi) is the reduced mass of two 
opposite-spin particles. Due to the fermionic antisymme¬ 
try, the regular parts are not independent functions: 

Ai 3 = A 24 = = —A 23 = A (2) 


Schrbdinger’s equation at zero eigenenergy A = 0, writ¬ 
ten in the language of distributions, is then 

= -4(r2 - R-13, i’4 - R-13)^(i’i - rs) 

- A{r2 - Ri 4, ra - Ri4)(5(ri - r4) 

- A{ri - R 23 , r4 - R23)i5(r2 - ra) 

-f A(ri - R24,r3 - R24)<5(r2 - r4) (3) 


with the kinetic energy Hamiltonian 


^ = E 


A 

-Ar 

2 m„ 


(4) 


and S(r) is the Dirac distribution in three dimensions, 
stemming from the identity Ar(l/r) = —47rS(r). 

We now go to momentum space and we take the 
Fourier transform of Schrodinger’s equation. In the left- 
hand side, each Laplace operator gives rise to a factor 
— where k„ is the wave-vector of particle number n. 
In the right-hand side, one obtains for example for the 


first term: 

4 


n<i 

n—1 


'3 


■'■"A(r2-Ri3, r4-Ri3)5(ri-r3) 


= (27r)3A(k2,k4)<5(^k„) (5) 


where the tilde indicates the Fourier transform. Intro¬ 
ducing the function D = (27r)^.A, we obtain the four- 
body momentum space ansatz generalizing to the 2-|-2 
fermionic problem the one of the 3-1-1 fermionic problem 


^tni(kl,k2,k3,k4) = 4 \ / [-P(k2,k4) 

E frkl 

- D(k 2 ,k 3 ) - D(ki,k 4 )-I-D(ki,k 3 )] ( 6 ) 


The ansatz obeys fermionic antisymmetry and 
Schrbdinger’s equation, not yet the contact condition o, 
that it suffices to implement for (i,j) = (1,3). One thus 
takes the inverse Fourier transform of at (ri, r 2 , r 3 , r 4 ), 
with the parametrization: 


ri = Ri3 + 

rs = Ri3 - 


m 3 

mi + m 3 
mi 

mi + m 3 


ri3 

ri3 


(7) 

( 8 ) 


Only the contribution '024 of I?(k 2 ,k 4 ) to ip diverges for 
ri 3 —>■ 0; in that inverse Fourier transform, we then take 
Ki 3 = ki -l-ks, ki 3 = /4i3(ki/mi -k3/m3) and k2,k4 as 
integration variables (clearly /ii 3 = nn), so that ki -ri -|- 

fi^kl tv^kl k?kl^ 

k 3 • r 3 = Ki 3 • Ri 3 -k ki 3 • ri 3 and -—- + -= -—- + 

2mi 2 m 3 2 fii 3 

— -311—integration over K 13 is straightforward, due 

2(toi -I- m 3 ) 

to the momentum conservation, and integration over ki 3 
also can be done using 


u(r) 

One obtains 


d^ki3 

(27r)3 kf^ + 47 rr 


(9) 


'024(ri,r2,r3,i'4) = J 


(Pk2cp'k4 2/413 


u{ri3) 


X e 

with qi 3 > 0 such that 

^^<li3 _ ^^(k 2 + k4)2 


(27r)® h? 

j [k 2 ■ (ra-Ri 3 )-l-k 4 ■ (r4-R 13 )] ^ 




h^kf 


2 /ii 3 2 (mi-|-m 3 ) 2 m 2 2m4 


( 10 ) 


( 11 ) 


Taking ri 3 —> 0 in '024 is then elementary. In the con¬ 
tribution to 0 of Il(k 2 ,k 3 ), Z)(ki,k 4 ) and D(ki,k 3 ), 



















4 


noted as ' 05 ^ 24 , one can directly take ri 3 = 0. Thanks 
to momentum conservation one can replace ki + k 3 by 
— (k 2 + k 4 ) within the position-dependent phase factor, 
so that the positions r 2 — R 13 and r 4 — R 13 appear as in 
Eq. (fTOl) : 


'0#24(ri = Ri3,r2,r3 = Ri3,r4) = 

[ d^fc2d^fc4 i[k2-(r2-Ri3)-Hc4-(r4-Ri3)] [ 

J (27r)9 J (27r)3 

[-^(k2,k3) - D(ki,k4) +i?(ki,k3)] (12) 

Finally the contact condition at the unitary limit, that 
is for 1 /a = 0 , leads to the following integral equation for 
D: 


3/2 


0 = 




(k 2 k4)^ 


kl 


2TTh? TO-|- -|- mi mi mi 


P 11/2 

^ i?(k2,k4) 






fd^kicfks 7 

(27r)3 


[i:)(k2,k3)-p:)(ki,k4)-i:)(ki,k3)] 


E 


2mr, 


(13) 


1313 

where the first term is simply Zl(k 2 ,k 4 ). Contrar- 

ily to the 3-1-1 fermionic case [l3,lQ, D is not subjected 
to any condition of exchange symmetry. 


III. TAKING ADVANTAGE OF SYMMETRIES 
A. Overview 

The unknown function Zl(k 2 ,k 4 ) in the integral equa¬ 
tion (US depends on six real variables. This is already 
a strong reduction, as compared to the twelve real vari¬ 
ables of the original four-body wavefunction, but still this 
makes a numerical solution challenging. 

Fortunately one can use rotational invariance as in sec¬ 
tion 1111^ the unknown function D can be considered for 
example as being the = 0 component of a spinor of 
angular momentum £. Then it is known how the var¬ 
ious 2£ \ components of the spinor transform under 

an arbitrary common rotation of k 2 and k 4 , in terms of 
rotation matrices having spherical harmonics as matrix 
elements, so that it suffices to know the value of the 2t+\ 
component of the spinor in the particular configuration 
where vector k 2 points along x axis in the positive direc¬ 
tion and k 4 lies in the xy upper half-plane y > 0 , at an 
angle 6*24 S [ 0 , 7 r] with respect to k 2 . As this particular 
configuration is characterized by the cosine of the angle 


024 and the two moduli p and p, the unknown function 
D(\i 2 , k 4 ) can be represented in terms of 2 £-|- 1 unknown 
functions fm] of these three real variables [T^ : 


e 

D{k2,k4) = ^ [F™"(e2 • e2,e4_L2 ■ 6^,624 ■ e^)]* 

mz ——£ 

X /m^(Pp4,U24) (14) 

In this expression we have introduced the unit vectors 


62 


64^2 

624 


k 2 

k2 

1 

^^24 




— ■a24e2 


k2 A k4 
|k2 A k4| 


(15) 

(16) 
(17) 


Here 024 G [0, tt] is the angle between k 2 and k 4 , and the 
notations 


M24 = cos 024 and V24 = sin924 (18) 

will be used throughout the paper. It is apparent that 
64 x 2 is obtained by projecting 64 = k 4 /fc 4 orthogonally 
to 62 and by renormalizing the result to unity. Then 
(e 2 , e 4 _L 2 , 624 ) forms a direct orthonormal basis. In that 
basis, an arbitrary (unit) vector n has uniquely defined 
spherical coordinates, that is polar angle 0n € [0, tt] with 
respect to axis 624 and azimuthal angle G [Oj 27r[ in 
the 02 — e 4 _L 2 plane with respect to axis 62 . Then 


y/"" (02 • n, 04 _L 2 • n, 024 • n) = Y™"(0„, (/„) (19) 


where the right-hand side is the standard notation for 
the spherical harmonics . Integral equations can then 


be obtained for the fm], see section IIII Bl 

For an infinite s-wave scattering length the Wigner- 
Bethe-Peierls contact conditions are scale invariant. 
As the integral equation (fT51) was further specialized to 
the zero energy case, its solution can be taken as scale 
invariant, which allows one to eliminate one more variable 


fml{k2,k4,U24) = (k^ + kj) (ch0;)®+^/^ 

X e*'"'^="/24.Wja;,a24) (20) 

with 

a; = In ( 21 ) 

k2 

The first factor contains the scaling exponent of the so¬ 
lution, which involves the unknown quantity s. By in¬ 
serting the ansatz (EUl) into the linear integral equations 
of section PlI Bl one obtains linear integral equations for 
the unknown functions 4>ml(a;,M), represented by a ma¬ 
trix M^^\s) that depends parametrically on s, see sec¬ 
tion [TTrCl requiring that the functions u) are not 
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identically zero, one gets an implicit equation for s, in 
the form 


detM(^)(s) = 0 ( 22 ) 

The way the first factor in Eq. o is parametrized by 
the quantity s ensures compatibility with the notation 
used by Efimov for the three-body problem Q. In the 
three-body problem, the Efimov effect takes place if and 
only if one of the scaling exponents s is purely imaginary, 
and the geometric trimer energy spectrum has then a ra¬ 
tio exp(—27r/|s|). In the four-body problem, with our 
definition of s, the four-body Efimov effect occurs if and 
only if there is a purely imaginary s solving Eq. (I22[) . in 
which case there exists a geometric sequence of tetramer 
eigenenergies with a ratio exp(—27r/|s|). A justification 
is given in [l3l| . The second factor in the ansatz (EOl) en¬ 
sures that the matrix M{s) is hermitian for purely imagi¬ 
nary s, with bounded diagonal matrix elements, which is 
both mathematically and numerically advantageous; as 
compared to [l^ it contains a additional term s in the 
exponent, that for purely imaginary s suppresses phase 
oscillations in the matrix elements of M{s) at large |a;| 
[50l| . The third factor in Eq. (1201) is a phase factor tak¬ 
ing into account the fact that exchanging k 2 and k 4 in 
Eq. (1141) transforms the spherical harmonics into 
(— with the same values of the variables 
[iSj; it ensures that the matrix M{s) transforms in the 
simplest way under the exchange of m-\- and mj,, which 
must leave our 2 - 1-2 problem invariant. 

A last reduction of the problem can be obtained from 
parity invariance. It turns out that, under the transfor¬ 
mation (k 2 ,k 4 ) —^ (—k 2 ,—k 4 ), the term of index rriz in 
the sum m acquires a factor (—1)™* [I^. This shows 
that the odd-parity functions (that is with niz odd) 
are decoupled from the even-parity functions <I>ml (that is 
with rriz even) in the integral equations, and that M^^\s) 
has zero matrix elements between the odd and the even 
channels. 


B. Rotational invariance 


To obtain the integral equations for the unknown func¬ 
tions fml in Eq. da we use a variational formulation: 
The integral equation (1131) is equivalent to 


dDHk2MAD,D*]=0 (23) 

where D and its complex conjugate D* are taken as inde¬ 
pendent variables, do* is the functional derivative with 
respect to D* and the functional £ is given by 


(24) 


with the diagonal part 


fdiag = J d^fc2d^/c4A’*(k2,k4)ZI(k2,k4) 


3/2 




(k 2 


2TTh? mi -|- m| 
and the generic off-diagonal part 


mi 


kl 


1/2 


(25) 


£24,ij — 


(P k2(P k4(P ki(P ks 


£>*(k 2 ,k 4 )Z/(ki,kj) 


(2^)3 

(5(ki -|- k 2 -l- kg -|- k 4 ) 

(fc| -f kl) 






(26) 


Then one inserts the ansatz m into these functionals. 
Assuming that one is able to integrate over all variables 
other than ^ 2 ,^ 4,^24 and ki,kj,6ij, one obtains a func¬ 
tional of the fml and fml*, which it remains to differenti- 
(i)* 

ate with respect to fmi to obtain the integral equations 
for the flnl ■ 

Integration is simplified by the following remark: The 

(i) 

final integral equations and their solutions fmi cannot 
depend on the specific vector introduced in Eq. da. 
One can then replace by an arbitrary unit vector n 
in the ansatz da, and one can average the resulting 
functional £ over n uniformly on the unit sphere, for fixed 
fmi- The result of this average is particularly simple 
when the orthonormal basis of Eqs. (I15I16I17I) reduces to 
the usual Cartesian basis: 


(02, e 4 _L 2 ,624) — (*^a: ; 7 *^2) 

Then [U 


(27) 


(E™" ( 62 -n, 64 X 2 -n, 624 -n))!),™" (exn, e|x*-n, 6 ^- 11 )] *)„ 
= ^((£,m,|i?to)|£,m;))* (28) 

where (•. •)n indicates the average over the direction of n 
and the quantum operator represents (in the usual 
spin-£ irreducible representation, with vectors |£, m^) of 
angular momentum rrizh along z) the unique real space 
rotation that maps the Cartesian basis onto the 

basis (e 2 ,e|X 2 ,*^ 2 /). 

(Oj , (29) 

After average over n, and integration over ki and kg in 
£ 24 , 1 ] I it remains an integral over k 2 and k 4 , with an 
integrand invariant by common rotation of k 2 and k 4 . 
To evaluate that integrand, one can then indeed assume 
that k 2 is along x (in the positive direction) and that k 4 
lies in the plane xy in the upper half y > 0: 

k2 = k2ex (30) 

k4 = A:4(cos 6(2462; -I- sin6(246y) with 624 G [0,7r] (31) 


£ — £diag + £24,23 + ^24,14 — ^24,13 
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so that 


( 02 , 64^2 ; ^ 24 ) — i^^XT^yi^z) 

and one can use Eqs. (I28I29I) . Then one pulls out a factor 
47 r (resulting from integration over the solid angle of k 2 ) 
compensated by the dvr denominator in Eq. (1281) , and an 
uncompensated factor 27r (resulting from the integration 
over the azimuthal angle of k 4 for the spherical coordi¬ 
nates of polar axis k 2 /fc 2 = e^; for k 4 ), and one is left 
with an integration over the moduli ^2 and fc 4 and over 
the angle 024 . 

For the functional Sdiag, this gives a simple result: 
Since i = 2 and j = 4, the matrix is the iden¬ 

tity matrix, i?hi) reduces to the identity operator; also 
there is no ki or ks integration. One obtains 


with U 23 = cos 023 and ti 23 = sin 023 as in Eq. (1181) . This 
gives 


^ poo pi pZTZ 

£ 2423 = Stt / dk 2 dk 3 dkik 2 klk\ / du23du24 dcj) 
(f, m'^)f^*{k2, k^, M24)/m- (* 2 ,^ 3 , U 23 ) 


(2^)3 


n?{k^+k‘^) . h?(k^+kl 

2 m'f- ' 2m^ 


(37) 


where fci is given by Eq. (1361) and we used the fact that 
has real matrix elements in the standard \l^ m^) basis. 


2. Case (i,j) ^ (1,4) 


£diag — 


E 


27r / dk2k2dk/ik\ 


I dU24\frr 


,(k2,k4,U24)\‘' 


f kl + kl + 2k2kiU24 ^ kj ^ fc| \ 
27:^? \ m| -I- mi m-f mi J 


with the same notation as in Eq. dni). For the generic 
off-diagonal part this leads to 


^ 24,ij — 


t «c 

I - D *^0 


dk2k2dkik^ / du24 


rriz ,m'— — 


c 

d 


d^kidd' fea 
(27r)3 

{{i, \e, m'S) f^l*{k2,k4, U24)f^l (h, kj,Uij) 


<^(ki -I- k2 -I- ka -I- k4) 
h4{k1+kl) 

' 2 m^ 


The way to proceed with the integration over the direc¬ 
tions of ki and ka depends on the indices i and j. 


For (i,j) = (1,4), one integrates over ka using the 
Dirac distribution, that imposes ka = — (ki -I- k 2 -I- k 4 ) 
and one integrates over ki using spherical coordinates in 
a rotated basis 

k4 

{ex,eY,ez) = ( 0^,64 A 0 ^, 04 ) with 64 = —. (38) 

ki 

The direction ez of k 4 is taken as the polar axis, so that 
the polar angle is 0 i 4 ; ex is taken as the azimuthal axis, 
with the azimuthal angle called (p, see Fig. [TJb). Then 
the real space rotation in Eq. (|2^ is 

= nz{<p-^)n4e24-6ii) 

= 7^.(024)7^.(0 - |)7^,(-0l4) (39) 

and the corresponding operator has matrix elements 

(i,m,\R^^‘^^\i,m',) = 6 "™^®=^" 

X (40) 


1. Case (i,j) = (2,1) 

For (i,j) = (2,3), one trivially integrates over ki using 
the Dirac distribution, that imposes ki = —(k 2 -l-k 3 -(-k 4 ), 
and one integrates over ks using spherical coordinates of 
polar axis e^ and of azimuthal axis e^; the azimuthal 
angle is called </>, and the polar angle is called 023 since 
it is the angle between k 2 and ka [see Fig. (Ha)]. Then 
Tliij) 

in Eq. (E^ is the rotation of axis x and of angle </>: 

7^(23) = and (35) 

where L^, is the angular momentum operator along x. 
Also 

^4 = kd ^ -(- /tg -)- k ^ -)- 2k2k3U23 -t- 2^:2^:41124 

-I- 2k3ki(u23U24 + V 23 V 24 cos (/>) (36) 


Using Bj; = [77.^^"^)] ^ 01 = U24ez -I- 1 : 246 ^ and 0 i = 
Ui4ez + vu(cos(j)ex -I- sin^ 0 y) we get 

fcg = fci -|- ^2 -f ^4 -I- 2 ki k2 [ltl 4 M 24 + V14V24 cos((f> — 7r/2)] 
-|-2/01/041114-|- 2 A: 2 A: 4 U 24 . ( 41 ) 


This gives: 


^ poo pi pZTT 

£24 14 = yy ‘2.17 / d/oic?/o 2 C?/o 4 /oi/o|/e 4 / duiidu24 / dp 
-^0 -^-1 -^0 


(27r)3 


h?{k'(+k'^) 

2m^ 


h?{k'^+k1) 


2m^ 


X /W*(fe 2 , ki, U2i)fy} (ki,ki, 1114 ) (42) 


where /03 is given by Eq. m- 
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FIG. 1: (Color online) Positions and parametrisations of the 
wavevectors appearing in the angular integration in the func¬ 
tionals S 24 ,ij ■ The vectors k 2 and k 4 are given by Eqs. (I30I31II . 
(a) For {i,j) = (2,3), ki = —(k 2 -I- ks -|- k 4 ) and one inte¬ 
grates over ka using spherical coordinates of polar axis x and 
azimuthal axis y. (b) For {i,j) = (1,4), ka = — (ki-|-ka-I-k 4 ) 
and one integrates over ki using the polar axis Z (direction of 
k 4 ) and the azimuthal axis X (direction of ej) as defined by 
Eq. (I38j, leading to the polar angle On and the azimuthal an¬ 
gle ((> (< 0 in the figure). The dashed line gives the direction 
of the component of ei = ki/fci in the XY plane, (c) For 
(*>j) = (li3), one integrates over the rotation TZ, such that 
ki and ka are given by the action of TZ on vectors kf^ and kf* 
in the xy plane as in Eqs. (I43I44II . using the parametrisation 
in terms of Euler angles associated to the convenient axes X, 
Y and Z of Eq. (|48)) . 


S. Case {i,j) = (1,3) 


For (i,j) = (1,3), we find it convenient to replace the 
integration over ki and ka by an integration over the 
moduli ki and fca, over the angle 9 13 G [0,7r] and over a 
rotation matrix TZ uniformly distributed over the rotation 
group SO(3), the vectors ki and ka being given by the 
action of TZ on vectors fixed in the xy plane: 


ki = 7^kf^ with kf’' = fcie,, (43) 

k 3 = 77 .k 3 ’" with = ks^uise^ + visSy) (44) 


Then TZ is precisely the rotation matrix of Eq. (|2^ 
and 


^^ 24,13 = ^ 


rriz ,m'— — 



n dknkl 


dui3du24 / dTZ 
-1 dsO(3) 


((£, m^\R\T, m(,))* f^l*{k 2 , ^ 4 , U 24 ) (fci, fca, M 13 ) 
<5(k2 +k4+7^(kf + kf")) 


2 m-|- ' 2 m^ 


(45) 


where the factor 2 originates from ( 47 r x 27r)^/[47r(27r)^], 
R is the operator representing TZ and dTZ is the invariant 
measure over the group SO (3) normalized to unity (see 
§8.2 of reference [^) [^. To integrate over TZ, we use 
the Euler parametrisation as in Eq. (7.1-12) of reference 


TZ = TZz{a)TZY{mz{l) (46) 


where the Euler angles a and 7 run over an interval of 
length 27r and the Euler angle (3 runs over [0, tt] , so that 
the invariant measure is (see § 8.2 of reference |^ ) 


da sin pd^d'y 

=- 8 P- 


(47) 


Due to the occurrence of k 2 -I- k 4 in the argument of 
the Dirac distribution in Eq. (l45)l . the convenient direct 
orthonormal basis defining the rotation axes X, Y and 
Z is now [see Fig. [Il[c)] 


/ N / . k2 -b k4 k2 -f k4 

= (48) 

Then the Dirac distribution can be written as 15311 


5(k2-l-k4-|-77.(k4''-|-k3’')) = (5(sin7)(5(sin(/3o-/3cos7)) 
^ ' 5 (|k 2 -b k 4 | cos(/3o - /3cos 7 ) -b |k^^ + k|^|) 

I sin/3o| jk^’' -b kf^'j |k 2 -b k 4 | 

where we have introduced the oriented angle /3o between 
kf’' -b kf’' and k 2 -b k 4 such that [see Fig. [Tl[c)] 

k®’^ -b ka’' = [kf’' -b ka’^K— sin/3o ex + cos/3o ez). (50) 

There is no dependence on a in the right-hand side of 
Eq. (|49|): in the argument of (5, one can write k 2 -b k 4 as 
'TZz{oL){k 2 + k 4 ) and, due to the rotational invariance of 
the three-dimensional Dirac distribution, one can factor 
out and remove the rotation TZz{a). The integration over 
a in Eq. (1451) then pulls out in the matrix element of R 
the orthogonal projector on the state of total angular 
momentum I and of vanishing angular momentum along 
Z-. 


/*27r 

/ = 27r|^, = 0 )(^, = 0 | 

JQ 


(51) 


























In the integral over 7, for example over the interval 
[—7r/2, 37 r/ 2 ], only the points 7 = 0 and 7 = tt con¬ 
tribute. The contribution of 7 = tt can be deduced 
from the one of 7 = 0 by changing /3 into —/3, due to 
= 7?.v(—/3) and to the invariance of 
k®’' + kfand | mz = 0) by rotation of axis Z. In the 
integral over /3 £ [0,7r], the 7 = tt contribution can then 
be taken into account by extending the integration of the 
7 = 0 contribution to /3 £ [— tt, 0]: one can take 7 = 0 in 
Eq. and one faces 


f fi/3| sin/3|e*^™^(5(sin(;3o-/3))'J(|k2+k4| cos{l3o-P) 
J —TT 

+ |kf + k«’^|) = <5(|k2 + k4| - |k«=^ + kf I) 

X f d/3| sin^ S{f3 — Pq — tt — 27rn) 

n&I. 

= 5(|k2 + k4| - |kf -f k«’'|)| sin^o|(-l)'"^e''''>™'^ (52) 

Due to the 27r periodicity of the integrand we have shifted 
the domain of integration so as to only keep for example 
the term n = 0 of the Dirac comb. Finally, using /3o = 
'r24 — ti 3 , where T 24 is the angle £ [ 0 , 7 r] between k 2 and 
k 2 -I- k 4 and T 13 is the angle £ [0, tt] between k®’' and 
kf’^ -I- kf’^ so that (up to a phase factor) 

|£, mz = 0) = TO,;, = 0), (53) 

and using the property that 

{£, rrix = 0|£, to(,) = 0 if ^ + to), is odd, (54) 
allowing one to replace (— 1 )'"* with (— 1 )^, we obtain 


^^24.13 = - / TT dknk^ / dui3dU24 

^ tf-o 27r 7o 


=0){i,mx =0 


k 2 -b k 4 

h^{k'^-\-k^) h^{k'^-\-k^) 

2m^ 2m^ 

|kfi’'-bkf='| 


X (5 (I k2 -bk41 -1 kf’^-bkfI) / (A:2 , ^4, U 24 ) fci, ^3, Mi3 ) 

(55) 


knowing that \^, rux = 0 ) has real components in the basis 
\£, rriz) up to a global phase, and that k 2 and k 4 are given 
by Eqs. (IMl^ and kf’^ and kf’^ by Eqs. (14311411 . 


C. Scale invariance 

To take advantage of the scale invariance of the zero- 
energy solution, one uses the ansatz ( 1 ^ with s £ zR, as 
physically explained in section IlII A[ and one inserts it in 
the various terms (I33I37I42I55I) of the functional (1^^ . In 
Eq. (1331) one performs in the integral over ^4 the change 
of variable ki = e^k 2 , where x ranges from —00 to -boo, 


and one sets U 24 
mass ratio 


u for conciseness, also introducing the 


rriT 
a = — 




(56) 


One pulls out a constant factor that will be given and 
discussed later, to obtain 


-'diag 


= T dx du 


1 L(i + «)^ 
1/2 




-b 


2 (a -b 1 ) cha;_ 


chx/ 
|$W(x,u)p (57) 


In Eq. (1571) one performs the change of variable ^4 = e^k 2 
and ^3 = ^2 in the integrals over ^4 and k^, also setting 

024 = 0, U 24 = u, V 24 = V and 023 = 0'■, U 23 = u', V 23 = v' 
for concision. One then pulls out the same factor T to 
obtain 


^ 24,23 = d- / dxdx / dudu 


rriz 


c,x-\-x' 


^chx' 

chx 


s/2 


1/4 




\4:chxch.x' 


/o 


(27r)2 'D24,23i4>',x,u]x',u';a) 


(58) 


In the denominator, we have introduced the notation 

^2/1.2 I i,2\ fr2{j.2 I 7,2 \ 


^24,23 — 


hr{ki+kl) hr{kl-\-ki) 

2m^ 2m^ 

Mt4 


(59) 


where ki is given by Eq. dMl) so that 


7^24,23(<(>; X, u; x', u'; a) = ch(a: — x') 

-\ - - —(e““““ -b e~^ u -b e~^u' -b uu' + vv' coscj)). 

1 + a 

(60) 

In Eq. (1421) one performs the change of variables k 4 = 
e^k 2 and ki = k 2 (so that ^ 4/^1 = ) in the inte¬ 

grals over k 4 and fci, and the change of variable (j) = ^+ 4 >' 
in the integral over <f) [ 5 ^, also setting 024 = 0 , U 24 = u, 
V 24 = V and 014 = 0', Ui 4 = u', V 14 = v'. Again pulling 
out the factor /F one gets 


^24,14 = F ^ 

ruz 


dxdx' / dudu' 

1 J-i 


e ch a; 


s/2 


o — (x+x') 


1/4 


ix',u') 


f 


4 ch X ch x' 


( 2^)2 


T’ 24 ,i 4 (</i’'; X, u; x', u'; a) 


( 61 ) 
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In the denominator we have introduced the notation 

^^2 / 1.2 I 7 2 \ Pi‘2 , 


V' 


24,14 = 


, h^(k‘+ki) 
2m-^ 2m^ 

h^kik2 


(62) 


with fca given by Eq. gB) so that 


2 ^ 24 , X, u] x', u'] a) = ch(x — x') 

H- ——u + e^u'+ uu'+ vv'coscj)'). (63) 

1 + a 

Finally, in Eq. (1551) . one performs the change of vari¬ 
ables ki = e^k 2 and k^ = fci in the integrals over 
ki and k^, also setting 024 = M 24 = u, r 24 = r and 

013 = 0', Ui 3 = u', ri 3 = r'. The integration over ki is 
straightforward due to the occurrence of a Dirac distri¬ 
bution in Eq. (l55l) . Due to the phase factor in the ansatz 
(EOl), there naturally appear the angles 7 = t — 0/2 and 
7 ' = t' — O' 12. Since r is the angle between k 2 and 
k 2 -I- k 4 [see Fig.[T][c)], one has according to Eqs. (I30I31I) 
and using the usual representation of vectors in the xy 
plane by complex numbers: 


e 


27 _ 


1 _l_ 

k + e e ^-i 9/2 
jl + e“e*®| 


Jx+ie)/2 I -{x+i6)/2 

- _ — _ (641 

|g(a;-|-je)/2 _|_ g-(a:-|-i6l)/2| 


As 0 G [ 0 , 7 r], the real part cosy of this expression is 
non negative so that one can choose 7 G [—f, §]■ Then 
forming the ratio of the imaginary part to the real part 
of the same expression gives the value of tan 7 and 

l + u) 

(65) 

One has the same expressions for 7 ', replacing the vari¬ 
ables X, 0 and u by x', 0', m'. This leads to 


7 = atan 


th 


(f) 


tan I - 


with tan ( — 


^24,13 — -A 


TTlz 


/ dxdx' / dudu' 
'r 0-1 


{u' + chx') chx' 


-| s /2 


{u -b ch x) ch X 


_ 2 _ 

47 r[('u + chx){u' + chx') chxchx']^/"^ 

{u + chx) + {W + chxO 

( 66 ) 


however that always the same diverging integral is pulled 
out, even if one singles out a wavenumber other than k 2 
(performing for example the change of variables /c 2 = 
e“^A :4 and fci = e““ ki in the integrals over k 2 and ki 
in Eq. (1421) 1. This is certainly due to the scale invari¬ 
ance of dk 2 /k 2 = d(hik 2 ). Alternatively, one can write 
the integral equation for the fmi deduced from the func¬ 
tional derivatives of Eqs. ()33I37I42I55|) of the functional 
Eq. with respect to fm *; at this stage, one has only 
used rotational invariance. Then, one inserts the scale 
invariant ansatz (EOl), and one obtains exactly the same 

it] 

integral equations for the as those derived from the 
functional derivatives of Eqs. (I57I58I61I66|) with respect 
to 


D. Parity invariance 


The term of index in the ansatz (1141) is simply mul¬ 
tiplied by (— 1 )™^ under the action of parity (k 2 ,k 4 ) — 
(—k 2 ,—k 4 ) [l 8 |. This means that the odd-niz compo- 
nents of are decoupled from the even-mz compo- 
nents of in the integral equation. This property can 
also be obtained by an explicit calculation: first, for 
and of different parities, the coupling amplitude be¬ 
tween and \£,m'^) must vanish in Eq. (1551) : this 

can be seen from Eq. (l54ll . Second, it also vanishes in 
Eqs. (j58l61|) after integration over </) or </>': obeys the 

selection rule Aruz = ± 1 , and in an expansion of 
in powers of </, only even powers of (j) and survive 
due to the parity of the denominator T). In what follows, 
at a given angular momentum £, we shall distinguish the 
manifold of parity (—l)^"*"^, where £24,13 and the contri¬ 
bution of D(ki, k 3 ) in Eq. (IT^ are zero, and the manifold 
of parity (—1)^ where they are a priori non zero. Note 
that, in the particular case £ = 0, there exists only the 
manifold of parity (— 1 )^- 


IV. FINAL FORM AND ASYMPTOTIC 
ANALYSIS 


A. Explicit form of the integral equation 


In all the results (j57l58l61l66l) there appears a factor 


Jo k 2 ■ 


(67) 


This factor contains a diverging integral, making these 
last calculations not entirely rigorous. We have checked 


By taking the functional derivative of £ of Eq. dm) 
with respect to , using the forms (I57I58I61I66I) of the 
various terms and not forgetting the minus sign in front 
of the last contribution in Eq. (IMl) . we obtain the form 
of the integral equation dm maximally reduced by use 
of the rotational symmetry and of the scale invariance: 
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0 = 


a 




1/2 


_(1 + q;)^V chx/ 2(Q; + l)cha;_ 
with the following expression for the matrix kernel 


f dx' f du 
Jk. J—1 


( 68 ) 




ch x' 


,\s /2 


oX+X 


1/4 


r27r 


d<j) 




e^'chx J \4cha;chx'y Jq (27r)2 ch(a; — a;') + y:^[(u + e“®)(M'+ e““') + hi;' cos 


' ch: 


1/4 




e ch X , 

(- 1 )^ 


4chxchx'y Jq (27r)2 ch(x — x') + y^[(m + e^)('«'+ e“') + wx'cos^] 

’ {u' + chx') chx'\ ma: = 0)(£, = 0|^, 7712 ) 6 “*™*'’' 


47r[(u ■ 


+ chx) ( 71 ' + chx') chx chx'] \ (u + chx)chx / ^ e-°^'+ae*’' ^ I I | 


(69) 


Here, the scaling exponent s is purely imaginary, so that a 
four-body Efimov takes place in our 2-1-2 fermionic prob¬ 
lem if Eq. (1551) has a non-identically zero solution for 
some non-zero s. We recall that the angle 0 G [0,7r] is 
such that u = cos 9 and 7 ; = (1 — = sin0, and that 

the angle 7 is given by Eq. (1651) : the same relations hold 
among the primed variables. 

The first, second and third contributions in Eq. (16911 
originate respectively from the terms Il(k 2 ,k 3 ), 
Zl(ki,k 4 ) and Zl(ki,k 3 ) in the unreduced integral 
equation m-, the diagonal term in Eq. (l 68 )) emanates 
from the diagonal term of that equation. The integrals 
over (j) can be evaluated after insertion of a closure 
relation in the eigenbasis of (s^ . Importantly, the 
third contribution in Eq. (1691) vanishes when I niz or 
i + ni'z are odd, i. e. in the parity channel (— 1 )^“''^, as 
shown by the property (1541) and as already pointed out 
in section UlIDl 

It is interesting to note the decoupled form of the 
prefactors in each contribution of Eq. (1551) . of the form 
[/(x, u)Y/^+^/yf{x', 77 ')]-''/ 2 +i /4 ^ith a function / given 
by 6^/(2 chx), e“^/(2chx) and l/[(77-|-chx) chx], respec¬ 
tively. The fact that this function /(x, u) is not common 
to all contributions prevents one from suppressing the s 
dependence of the matrix kernel by a simple gauge 
transform on As expected, the s-dependence of the 
problem (16811 is non trivial. 

Our results (I68I69I) must obey the symmetry of the 
2-1-2 problem under the exchange of f and 4 ,. First, this 
exchange has the effect of changing the mass ratio a into 
its inverse 1/a, see Eq. (1551) . Second, the momenta k 2 
and k 4 in £)(k 2 , k 4 ) are exchanged, so that x of Eq. (l2T1) 
is changed into its opposite; this also reverts the direction 
of the quantization axis k 2 A k 4 along which the angular 
momentum rtiz is measured in Eq. dH): it changes rriz 


into —rriz according to the identity 


s--^-/''l£,m,) = (-l)'l^,-m,); 


(70) 


on the contrary, the non-oriented angle 6*24 € [ 0 , tt] be¬ 
tween k 2 and k 4 is unchanged, so that the variable u is 
unaffected. Hence, one must have 


a:', u'; a) = (-x, u\ -x', 77 '; a ^) 

(71) 

for all values of the argument and of the indices of the 
kernel. It is clear that Eq. (1591) indeed obeys the symme¬ 
try requirement (ED: the first and second contributions 
are interchanged, whereas the third one is invariant since 
7 is changed into — 7 , see Eq. dlS]). Note that our re¬ 
sults also respect the parity invariance, see section UlI D1 
and that the matrix kernel is hermitian as our variational 
derivation guarantees: 


ix,u;x',u';a) = ^^x',u';x,u;a) 


(72) 


B. Recovering the three-body problem from 
four-body asymptotics 

The right-hand side of the integral equation (155)) de¬ 
fines an operator M^^\s) acting on the spinor functions 
77 ). The spectrum of this operator is physically 
relevant, since a four-body Efimov effect takes place with 
an Efimov scaling exponent s € zR if and only if one of 
the eigenvalue Vt of is zero. As is a her¬ 

mitian operator, since s is here purely imaginary, its spec¬ 
trum is real and in general includes a discrete part and 
a continuous part. The discrete spectrum corresponds to 
localized, square integrable eigenfunctions; we are able 
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to determine it only numerically. The expected contribu¬ 
tion to the continuous spectrum corresponds to extended 
functions, that explore arbitrarily large values of \x\; as 
we now explain, it can be determined analytically from 
the asymptotic analysis of the kernel (1691) when x and 
x' tend to ±oo, by a generalisation of the discussion of 
reference M- There is also an unexpected contribution 
to the continuous spectrum, whose analysis is deferred to 
section II V Cl 

Sector X —>■ +00, x' —>■ +00: Clearly, the diagonal 
part of in Eq. (16811 tends exponentially rapidly 

to a finite and non-zero value, and the second and third 
contributions to the kernel in Eq. (1691) tend exponentially 
fast to zero. In the first contribution in Eq. (l69ll . the pref¬ 
actor tends exponentially to unity since e'“/(2chx) —)► 1, 
and in the denominator of the integrand, all the x or 
x'-dependent terms are exponentially suppressed, except 
the first one ch(x —x') since no hypothesis must be made 
on the magnitude of the difference x — x'. The eigenvalue 
problem then asymptotically reduces to 




Qf(2 -|- a) 

{l + a)\ 


1/2 „ i 

dx'Y^ 

m'=-e 



/o {2 ttY ch(x - xO + 


(x',m') 


(73) 


where the arrow in the exponent of 12 indicates that x and 
x' tend to positive infinity, and the phase factors 
and have been eliminated by a gauge transform 

on the spinor, <I’ml(x,'u) = (x, m). Then, 

one performs a spin rotation, by moving to the internal 
basis of the eigenstates \i,mx) of L^, in which jg 

diagonal: the components {x, u) are all decoupled. 
Eor a given rrix, with \mx\ < the trick is to extend 
(x, v) into a function of the real variable x and of 
the vector n on the two-dimensional unit sphere: 

EW(x,n) = $W(x,cosd)e™-^ (74) 

where 9 G [0,7r] and (/> G [0,27r] are the polar and az¬ 
imuthal angles of n in spherical coordinates, e.g. with re¬ 
spect to X and j/ axes, n = (cos0,sin0cos0,sin0sin;/)). 
In the phase factor in the numerator and in cos (j)' 

in the denominator, one can then replace by (j)’ — (j>: 
the integrand is a periodic function of of period 27r 
and its integral has the same value whatever the inter¬ 
val of length 27r over which </)' runs. Then, one recog¬ 
nizes the scalar product n • n' = uu' + vv' cos{(j) — cj)') 
where n' = (cos 0', sin 0'cos (()', sin 0'sin (/>'). The eigen¬ 
value problem is now 


1 -I- a 

F^l{x',n') 


dx' 


d'^n 


'|n|=i (27r)2 ch(x - x') -f 


nn 

l+a 


(75) 


The corresponding operator is invariant by translation 
along X and by rotation of n over the unit sphere. Its 
eigenfunctions Fml (x, n) can therefore be taken as plane 
waves of the variable x and spherical harmonics of the 
variables (0, </>), with the same quantum number m^, (this 
is imposed by the form (|7^ ) but with any integer quan¬ 
tum number L > \mx\ for the total angular momentum: 

FW(x,n)=e*'=“Yr-(0,0) (76) 

As usual for a rotationally invariant operator, the eigen¬ 
value does not depend on rux- It only depends on A, 
so it suffices to specialize to = 0, where Y]^{9,(j)) oc 
Pl(cos0), where Pl{X) is the Legendre polynomial of 
degree L. Then, one gets the continuous spectrum “to 
the right” (x,x' —>■ -l-oo): 

e {Ai(zA,a-i),V/c e R,VL > |m,|} (77) 

The function A^ of s G and of the mass ratio was 
introduced and analytically calculated in [H, , gener¬ 
alizing previous results [^, : 


. . (1 + 2 ^)'/' 
Al(s,^) = 


1 + ; 


du 


= cosv{P)+ — 




sini/(/3) 


dePL 


dx e ’^^Pl{u) 

27r chx -I- 

sin(s0) 


cost 


sinzz(/3) / sin(s7r) 

(78) 


where, in the second expression obtained after integration 
over X [ 3 ^ . 


v{P) = asin (79) 

is a mass angle. Eor all /3 > 0, it is found numerically for 
even L that the maximal value of Ki^{s,0) over s G 
is reached at s = 0, and the minimal value is reached 
for |s| —>• - 1-00 (where Al(s,/ 3) tends to cost 2 (/ 3 )). For 
odd L, the situation is found to be reversed: Al{s,/3) is 
minimal at s = 0 and maximal at infinity. To summarize, 
we expect that 

cosi^(/3) < Al(s,/ 3) < Ai(0,/3) Vs G iR,L even, 
Al( 0,/3) < Ai(s,/?) < cosz/(/3) Vs G iR,L odd (80) 

This allows to determine the borders of the continuous 
component of quantum number L in Eq. (1771) . see Fig. [7] 
A physical explanation for the emergence of the function 
Ai is postponed to the end of the section. 

Sector X —>■ —oc,x' —>■ —00: The calculation closely 
resembles the previous one, except that it is now the 
second contribution in the right-hand side of Eq. (15^ 
that survives. This was expected from the symmetry 
relation dZD. We arrive at the continuous spectrum “to 
the left” (x,x' —>• — 00 ): 

G {AL(zk,a),Vk G R,VL > \mx\} (81) 
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that differs from (113 by the occurrence of a (rather than 
1/a) in the argument of the function |57l |. The bor¬ 
ders of the L-components of that continuum are plotted 
in Fig. [3 for the first few values of L, using the numeri¬ 
cally checked property (l80l) . 

Parity considerations: At fixed £, the results (I77I81I) 
are expressed in terms of the quantum number 
whereas the original problem only distinguishes between 
an even parity manifold (m^ is even) and an odd parity 
manifold (m^ is odd). In practice, due to the property 
(IMl) . the continua (I77I81I) with L = 0 can be realized 
only in the manifold of parity (—1)^, at any considered 
total angular momentum I (obviously one must then take 
rux = 0). The other continua (with L > 1) can all be re¬ 
alized, in both odd and even manifolds, for all values of 

£ > 0 [H. 

Physical discussion: The function Ai(s,/3) appears in 
the unitary three-body problem of two fermionic parti¬ 
cles interacting with a single distinguishable particle, f3 
being the mass ratio of the majority-to-minority species. 
For s G jK. this function is given by the first form in 
Eq. (1751) : it can be analytically extended to real values 
of s using e.g. the second form in Eq. dUl) [13. The 
zero-energy solutions of this three-body problem have an 
Efimov scaling exponent s: the three-body wavefunction 
scales as R being the three-particle hyperradius, 

and the allowed values of s at total angular momentum 
L must solve 

Al(s,/3)=0. (82) 

This three-body system exhibits an Efimov effect if and 
only if this equation has a purely imaginary solution s € 
zR*. This occurs only at odd L, starting from a mass 
ratio (13 

/3 > ac(2;l) = 13.60696... (83) 

for L = 1, and at increasingly larger critical mass ratios 

for A = 3,5,... [H 113 . 

It is thus apparent that the asymptotic analysis of the 
2 + 2 fermionic problem brings up the three-body prob¬ 
lem. This is intuitive in position space: imagine that at 
fixed position r 4 7 ^ 0 of the fourth particle (of spin ),), the 
positions (ri)i<i <3 of the other particles (of spin ft/-) si¬ 
multaneously tend to zero; then the four-body wavefunc¬ 
tion ^(ri, r 2 , ra, r 4 ) must reproduce the behavior of the 
zero-energy scattering state of two f and one | particles, 
characterized by a mass ratio /3 = = a, in partic¬ 

ular it must exhibit the same scaling exponents s as the 
2 + 1 problem (see §5.3.6 in reference [S^); as these scal¬ 
ing exponents solve Eq. (1821) with /3 = a, this explains the 
occurrence of Al(s, a) in the spectrum (ISTll [^ . Even if 
£ = 0 for the four-body system, L can take any value, as 
the angular momentum can be distributed among parti¬ 
cle 4 and the first three particles. The equivalent in mo¬ 
mentum space of the considered limit is to have divergent 
(ki)i<i <3 at fixed k 4 , which, due to scale invariance, is 
equivalent to having k4 — > 0 at fixed (ki)i<i<3, that is 



FIG. 2: (Color online) Analytically obtained borders of the 
continuous spectrum of M{s) (s G zR’*") corresponding to the 
limits X ±00. It is a collection of components characterized 
by an angular momentum quantum number L of a three-body 
asymptotic problem (not to be confused with the total angular 
momentum i of the four-body states), leading to a continuous 
set of eigenvalues ranging from A_l(0, a) to cos[j/(q:)] for x —>■ 
—00, and ranging from Ai:,(0,1/a) to cos[j/(l/a)] for x —^ +00, 
where a is the mass ratio given by Eq. dSi), the An function 
is given by Eq. (1781) and the mass angle v is given by Eq. (17911 . 
We plot cos[:z(q)] (black solid line) and Ai(0, a) (color dashed 
lines, with L = 0,2,4 from top to bottom above the solid 
line, and L = 1, 3 from bottom to top below the solid line) as 
functions of a G [0, 0^(2; 1)] where ac(2; 1) is the critical mass 
ratio (1831) for the Efimov effect in the ttf three-body problem. 
Due to the a .f-J- 1/a symmetry of the 2 + 2 fermionic problem, 
one can restrict to a > 1 (that is to the right of the vertical 
dotted line); the borders of the x —^ —00 continuum can then 
be directly read on the figure, and the ones of the x +00 
continuum can be obtained by mentally folding back the a < 1 
part of the figure into the a > 1 part. 

X ^ —00 according to Eq. This is why /3 = a cor¬ 

responds to the spectrum (1811) . A similar reasoning with 
r 2 fixed with tending to zero leads to (3 = 1/a 

and X —> + 00 , as for the spectrum dZZl). 


C. A third, unexpected continuum 

The first two contributions in Eq. (1551) are innocuous: 
the denominator in their integrands cannot vanish, see 
Eqs. (I59l60p and Eqs. (162163^ . and, as we have seen, they 
have a short range in the (x,x') space. On the contrary, 
the third contribution in Eq. which is non-zero only 
in the (— 1 )^ parity sector, diverges when {x, u) —>■ ( 0 , — 1 ) 
or {x',u') —>■ (0,-1). This creates doubt about the 
bounded nature of the eigenvalues of M(s), s G zR, for 
that parity. We investigate this problem mathematically 
in the aopendixEland we conclude that M(s) is bounded. 

Physically, this divergence of the kernel leads to a quite 
interesting effect: the emergence of a third component of 
the continuous spectrum of M(s), different from the pre¬ 
viously discussed x —>■ ±00 continua. The intuitive idea 
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is that one can turn the eigenvalue problem r 2 <i> = M (s)$ 
into an integral equation with a bounded kernel through 
an appropriate change of variables, with the consequence 
that one of the new variables, that we shall call t, can 
tend to — oo, in which case the eigenvector $ takes a 
plane wave structure oc exp(ifct), fc G R, with a spec¬ 
trum: 

(84) 

This is an unexpected feature of the 2-1-2 fermionic prob¬ 
lem, absent in the 3-1-1 fermionic case [l^ . 

To obtain this result, we construct a local approxima¬ 
tion to the integral equation in the vicinity of {x, u) = 
( 0 , — 1 ), (x', u') = ( 0 ,- 1 ), keeping only the leading di¬ 
verging contributions. We use 


The resulting eigenvalue problem 




21/2 




(- 1 )^ dt'm 

2 ^ 1 ‘^-k ch{t - t') 


(92) 


admits Eq. (1841) as a continuous spectrum with eigen¬ 
functions (pit) that are for t —>■ —oo linear superpositions 
of and since / cht = tt/ ch(fc7r/2); we 

have checked numerically that it has no discrete eigen¬ 
value [6nl |. 

Physical interpretation: We collect Eqs. (I86I89I91L 
taking as a particular solution of Eq. (IMl) at t large and 
negative the function pit) = 1 , which corresponds to an 
asymptotic plane wave in t space with a vanishing wave 
vector, that is to fc = 0 in Eq. (IMl) [ 6 l|. Restrictin g fo r 
simplicity to a zero total angular momentum i = 0 [62j| . 
we then find that 


y = TT — 9 


(85) 


rather than u = cos 9 as integration variable, so that 
y,y' —^ 0 when u,u' —>■ —1. This pulls out a Jaco¬ 
bian sin y' that we absorb (in a way preserving the her- 
miticity of the problem) with a change of function. We 
also take into account the fact that the third, diverg¬ 
ing contribution in Eq. (1691) involves a projector onto the 
|£, rux = 0 ) state, and that the phase factors and 

^ [(y-i-chx) chx]“®/^ and [(w'-fchx') chx']'*/^ can 
be eliminated by a change of gauge. Hence the ansatz 


4*^1 [x, u) = 


(x^ -I- 2 /^) ®/^ 


,1/2 


= 0 )$(x,?/) 


( 86 ) 

where sinj/ was linearized, u-|-chx was quadratized. The 
resulting local eigenvalue problem is 



4 >(x, 2 /) 


(- 1 )^ [ dx'dy'jyy'f!^ 

2^/^TT Jjj [(x2+y2)(^/2+y/2)]l/4 

X „ (B^) 


where we have conveniently restricted the integration to 
the upper half {y > 0) of the disk D of radius po 1 
centered in (0,0). In polar coordinates 


{x,y) = {pcostp, psintp) ( 88 ) 

only {yy')^^^ depends on ip in the kernel. It depends on 
tp in a. factorised way so that $ is also factorised: 

$(x, 2 /) = p-^/^^ip)isini))^/^ (89) 

and, since dip' sin' 0 ' = 2 , 

1 \ (- 1 )^ 21/2 P 

n 2 I ./ 2 ^(P ) 


Si-7j)4-(rf = - 


P 


(90) 

The scale invariance of this kernel motivates the logarith¬ 
mic change of variable 

1/2 


t = In — and p{t) = 
PC 


-) 

PoJ 


$(p) (91) 


<I>®(x, u) oc 


{x,u)->-{0,-l) pS+3/2 


(93) 


A more inspiring writing is obtained in terms of the 
center-of-mass and relative wavevectors K 24 = k 2 -f k 4 
and k 24 = (k 2 — ak 4 )/(l -|- a) of particles 2 and 4: 


$^°^(x, u) oc 

K24jk2A^0 


K 24 ) 


(94) 


One has indeed p^ ~ 2(it -I- chx) and K 24 = 2fc|e^(M -I- 
chx), so that K 24 and p vanish in the same way when 
(x,m) —(0,-1); also the ratio ^^ 24/^24 tends to zero if 
and only if m -|- chx —>■ 0 [6^. Restricting to a small 
neighbourhood of the singularity, K 24 < efc 24 , where 
e <C 1 , we can in the ansatz dsni) approximate the factor 
(chx )®“'"^/2 by one and, in the denominator, approximate 
kl + kl = 2 fc ^4 -H 2^k24 ■ K 24 + KI 4 by its leading 

order approximation 2 k‘^A to isolate the singular behavior 
of 0 (k 2 ,k 4 ): 

^-g(k 2 ,k 4 )cx^(^£j (95) 

The key idea is then to see how this translates into 
a singularity of the regular part A 13 of the four-body 
wavefunction that appears in the Wigner-Bethe-Peierls 
contact condition ©■ As we have seen below Eq. ®, 
Ai 3 = a is related to II(k 2 , k 4 ) by a Fourier transform; 
using (k 24 ,K 24 ) rather than (k 2 ,k 4 ) as integration vari¬ 
ables, and the fact that k 2 -r 2 -l-k 4 -r 4 = k24-r24-l-K24-R24, 
where r24 = r 2 —r 4 and R 24 = ( 7712 X 2 -(- 7714 X 4 )/(m 2 -I- 7714 ) 
are the relative and center-of-mass coordinates of the par¬ 
ticles 2 and 4, we obtain for the contribution to A of the 
singularity of D: 


‘x4sing(x2 Ri3i 1*4 R 13 ) CK 


J i^24 <£^24 

pjk24 r24 / m . \ ^*+3/2 
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Integrating over the solid angles for k 24 and K 24 , per¬ 
forming the change of variable 1^24 = 9 fc 24 ''' 24 /|R 24 ~R-i 3 | 
at fixed ^ 24 , changing the order of integration over ^24 
and q and finally integrating over ^24 |64| we obtain 


'^sing(^^2 1^13; ^4 1^13) ^ 


R24-Rl3r~"/^ 

^24 


^e|R24 —R'13|/'^24 


dq 


qs4-l/2 


[|5_l|-l/2_(g+l)-l/2] 


(97) 


It then becomes obvious that the singularity in II(k 2 , k 4 ) 
at k 2 -l-k 4 = 0 is linked to a l/7’24 divergence of the regu- 
lar part ^13 of the four-body wavefunction at ^24 = 0 [65| 
(^ . This was physically expected: v 4 i 3 (r 2 — Ri 3 ,r 4 — 
R 13 ) is essentially the wavefunction of particles 2 and 
4 knowing that particles 1 and 3 have converged to the 
same location in the s-wave; since 2 and 4 are in differ¬ 
ent spin states, they interact in the s-wave and must be 
sensitive to the 2 — 4 Wigner-Bethe-Peierls contact con¬ 
ditions, which implies a l/r 24 divergence when r 24 —> 0 . 
Such a divergence of ^13 was already pointed out in the 
scattering problem of two fj- dimers in reference [4l| , and 
in the general fermionic problem when > 2 

and Ni > 2 in reference (see footnote 20 of that 
reference) [^ . 

This interpretation of the matrix kernel singularity at 
(x,u) = ( 0 , — 1 ) has a simple, though illuminating impli¬ 
cation: The 1 /r 24 divergence in ^13 can take place only 
when the particles 2 and 4 converge to the same location 
in the relative partial wave of zero angular momentum, 
since f and 4 - particles resonantly interact only in the s- 
wave. In such a configuration the angular momentum i of 
the function ^13 (that is of the whole system) is carried 
out by the center of mass motion of particles 2 and 4 with 
respect to R 13 ; in this case there is a univocal link be¬ 
tween the angular momentum i and the parity, as for sin¬ 
gle particle systems, and the parity of A 13 must be (—I)^- 
This explains why the singularity at (a;, u) = (0, —I), and 
ultimately the third continuum (1841) . can appear only in 
that parity channel 


V. SEARCH FOR THE FOUR-BODY EFIMOV 
EFFECT 


In the 3 -|- I fermionic problem, the signature of a 
four-body Efimov effect was that an eigenvalue of the 
corresponding M(s = 0 ) operator crosses zero for some 
value of a below 0 ,^( 2 ; 1) ~ 13.6069, specifically for 
a = q;c( 3; I) ~ 13.384 [Q. The question here is to know 
whether or not such a crossing can occur for the 2 - 1-2 
fermionic problem, that is for M(s = 0 ) corresponding 
to Eqs. (I68I69I) . We answer this question by a numerical 
calculation of the eigenvalues of M(s = 0 ). 

Numerical implementation: In the (-1)^+^ par¬ 
ity sector, we truncate the x variable in a symmet¬ 
ric way, that is to [imin = -a:max, a^max], and we dis¬ 


cretize it with a uniform step dx according to the 
usual mid-point integration method. We use 6 rather 
than u = cos 6 as integration variable, so that we use 
$(x, 0 ) = (sin 0 )^/^$(a;, u) rather than <I)(a;,M) as the 
unknown function; multiplying the eigenvalue problem 
124) = M{s = 0)4) by (sin0)^/^, we get an hermi- 
tian eigenvalue problem with the same eigenvalues, the 
same diagonal part and a kernel {x,0;x',9') = 

(sin 0 sin 0 ')^/^A^^ (a;,u;a;',u') with u = cos 0 , u' = 

COS0'. For better accuracy, we use the Gauss-Legendre 
integration scheme [ij for 6 ' as this angular variable is 
naturally bounded to [ 0 , 7 r] [^ . 

In the (—1)^ parity sector, there is the additional com¬ 
plication that the kernel diverges at the point {x,9) = 
( 0 , 7 r), making the previous {x,9) discretisation inefh- 
cient. At a small distance from this point, say less than 
po, the optimal set of variables is where ijj G [ 0 , 7 r] 

is defined by Eq. (1551) and t G is given by Eq. (ITO . 
since the resulting kernel is bounded after a convenient 
change of the unknown function (t,^) = p^iil{x,9), 
see Eq. dMl). So we resort to a mixed scheme: for 
p = [x^ + (tt — 0)^]^/^ > Po, we use the (x,9) set 
of variables, with x uniformly discretised with a step 
dx submultiple of po and 6 discretised according to a 
Gauss-Legendre scheme over the interval [0,0niax] where 
^'max = Tr for |a;| > po and 0niax = - (Po - for 

1 2 ^I < Po, the number of angular points being scaled lin¬ 
early with Sinax) for p < Po, we use the (t,'0) set of vari¬ 
ables, with t truncated to [tmin, 0 ] and discretized with 
a uniform step dt according to the mid-point integra¬ 
tion method, and '0 G [0, tt] discretised according to the 
Gauss-Legendre scheme [70l |. 

Results: The numerically obtained spectrum of M(s = 
0) for the 2-1-2 fermionic problem is plotted in Fig. [3] 
[left half for the (—1)^ parity channels, right half for the 
(—1)^+^ parity channels], for the first values 0 < ^ < 3 
of the four-body internal angular momentum quantum 
number as a function of the mass ratio a = —. It 
is a symmetric function under the exchange a O I/a so 
the figure is restricted to a > I; as the starting integral 
equation (|I3I1 assumes scale invariance, which is broken 
by the three-body Efimov effect beyond the threshold 
ac(2; I) = 13.6069 ..., the figure is also restricted to a < 
ac( 2 ;I). 

For the (—1)^ parity channels, the spectrum is entirely 
within the limits of the analytically predicted continu¬ 
ous spectrum, which are shown as dashed lines, except 
for even £ in a barely visible small triangle 0 close to 
U = 0.75 with I < a < 1.2: this means that the neces¬ 
sarily discrete numerical spectrum must tend to a con¬ 
tinuum when the truncations Xmax = — Xmin and tmin 
tend to -koo and —00 respectively. The lower border of 
the continuum corresponds for even i to the fc —> 0 limit 
in Eq. (15^ . that is zero, and for odd I to the smallest 
of the two quantities, I/-\/2 [this is the k —>• -|-oo limit 
of Eq. (15^ ] and Ai=i(0,a) [this is the minimal value of 
Eq. (1511) . see Fig. [5]. The upper border of the continuum 
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FIG. 3: (Color online) Eigenvalues of M(s = 0) as functions of the mass ratio a G [1,13.6], for various values of the angular 
momentum I and the parity [left half for the parity (—1)^, right half for the parity (—1)^”''^], obtained numerically after 
discretisation and truncation of the variables x and 6 in the zone p = [x^ + [it — 0)^]^^^ > po and of their log-polar versions 
t = ln(p/po) and ■(/) in the zone p < po (see text): imax = —^min = 12, dx = 1/5, dd ~ vr/lS, po = 2/5, tmin = —12, dt = 1/5, 
dtp = 7r/15, with the Gauss-Legendre integration scheme for the integrals over 9 and tp, and the mid-point integration formula 
for the integrals over x and t. For the (—l)^'*'^ parity channels the p < po zone is not useful and is not included in the numerics. 
The boundaries of the continuous spectrum of M(s) are given by black thick dashed curves: for the (—1)^”'"^ parity channels, 
this corresponds to the a; -I- ± oo continua of Eqs. (I77I81I) with all L > 1; for the (—1)^ parity channels, this corresponds to 
the a; —>■ -I- ± 00 continua of Eqs. (I77I81II with all L > 0 and to the {x,9) (0, tt) [that is t —>■ —oo] continuum of Eq. (1841) . 
Contrarily to the 3 + 1 case, no eigenvalue of M(s = 0) is found to cross zero for a < q:c( 2; 1) = 13.6069 ...: no four-body 
Efimov effect is found for the 2 + 2 fermionic problem. 


corresponds, whatever the parity of £, to Ai=o( 0 , l/a) 
[this is the maximal value of Eq. (1771) . see Fig. [5]. 

For the (—1)^“''^ parity channels, there are three dif¬ 
ferences. First, the continuum (IMl) cannot be realized, 
so that the lower border of the continuous spectrum of 
M{s = 0), now given by Al=i(0, a), reaches zero only at 
a = ac{2- 1). Second, the L = 0 continua in Eqs. (1771 ISTl) 
cannot be realized so that the upper border of the contin¬ 
uum of M(s = 0), given by Al= 2 ( 0 , l/a), is everywhere 
below 1 = limQ,_>_|_oo Al=2(0, l/a). Third, the contin¬ 
uum presents, in the a — plane for 1.53 ^ a, a large 
void internal area, corresponding to Ai= 2 ( 0 ,a) < 11 < 
Ai=i(0, l/a). Still, many numerically found eigenvalues 
lay in this internal area: these eigenvalues must corre¬ 
spond to the discrete spectrum of M(s = 0 ), with local¬ 
ized (square-integrable) eigenfunctions [t^. We checked 
this numerically by calculating the density of states of 
M{s = 0), in practice the histogram of its eigenval¬ 
ues, for increasing values of the numerical truncation 
a^max = — by doubling the values of Xmax and 

Xmin, the spacing m Tr/xmax between successive k in 
Eqs. (I77I81I) is approximately divided by two so that the 
density of states of the numerical quasi-continuum is ap¬ 
proximately multiplied by two, whereas the density of 
states of the discrete spectrum is essentially not (only 
exponentially weakly) affected as soon as imax is much 
larger than the localized eigenfunctions width in x space. 
This is what is observed in Fig. SI knowing that the lo¬ 


cations of the internal and external borders of the con¬ 
tinuum (corresponding to the dashed lines in Fig. [3|) are 
indicated by vertical dashed lines in Fig. S) 

Synthesis: Contrarily to the 3+1 case, no discrete eigen¬ 
value of Af(s = 0) (necessarily discrete because it would 
be below the lower border of the continuum) crosses zero 
for a < ac(2; 1) = 13.6069 ..., that is below the three- 
body Efimov effect threshold: no four-body Efimov ef¬ 
fect is found for the 2 + 2 fermionic problem [^. This 
conclusion is apparent in Fig. [3l obtained for all the in¬ 
ternal angular momentum quantum numbers 0 < £ < 3. 
It extends to all the angular momentum values that we 
were able to numerically explore, 4 < £ < 12, as we 
have shown with a dedicated careful spectral analysis 
almost perfectly at the three-body critical mass ratio, 
a = 13.6069, see Fig. [SJ 


VI. CONCLUSION 

We have studied in three dimensions a four-body 2 + 2 
fermionic system with resonant interactions and we have 
derived its momentum space integral equations at zero 
energy. By using rotational invariance and scale invari¬ 
ance, we have reduced them to a numerically tractable 
two-dimensional form (the unknown function depends on 
two variables only). With these equations we have nu¬ 
merically shown that no four-body Efimov effect occurs 
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FIG. 4: (Color online) Histogram of the eigenvalues Q of 
M(s = 0) for a mass ratio a = 10 and an angular momentum 
^ = 1 in the ( — 1)^'*'^ parity channels, obtained numerically af¬ 
ter discretisation and truncation of the x and 9 variables. The 
numerical grid is the same as in Fig.[3l except that much larger 
values of = —a^min are used, to reveal the emergence of 
the continuous part of the spectrum in the Xmax —> +oo limit: 
imax = 48 (red bars in the foreground) and Xmax = 96 (blue 
bars in the background). The black vertical dashed lines in¬ 
dicate the analytically predicted borders of the continuous 
spectrum (as in Fig. 0; in between the first two ones and 
in between the last two ones, it is indeed observed that the 
number of eigenvalues per bin is approximately multiplied by 
two when Xmax is doubled. On the contrary, the histogram 
is unaffected by the change of Xmax in the bins strictly in be¬ 
tween the second and third dashed lines, indicating that the 
corresponding eigenvalues belong to the discrete spectrum of 
M(s = 0), with localized eigenfunctions in x space. 


for the 2 - 1-2 fermionic system in angular momentum chan¬ 
nels 0 < £ < 12. The 3-1-1 fermionic system thus remains 
the only known one exhibiting a four-body Efimov effect 

0 - 

A detailed treatment of the second motivation for de¬ 
riving these integral equations, that is a calculation of 
the already measured [^, fourth cluster coefficient 64 
of the spin 1/2 unitary Fermi gas, is beyond the scope of 
this paper. Still we have numerically calculated an ed¬ 
ucated guess for 64 inspired from the analytical form of 
the third cluster coefficient 63 [H, : it is found that 

this guess does not reproduce the experimental value (see 
Appendix |B]) so that a dedicated work is needed and left 
for the future. 


Acknowledgments 

We acknowledge useful discussions with J. Levinsen 
and D. Gridnev. S.E. thanks the Japan Society for the 
Promotion of Science (JSPS) for support. The group of 
Y.C. is affiliated to IFRAF. 



1/Tmax’ 


FIG. 5: (Golor online) Numerically determined minimal 
eigenvalue flmin of M^^\s = 0) almost at the three-body 
critical mass ratio, a = 13.6069 ~ ac(2; 1), as a function 
of the numerical cut-off tmin = —Smax. For each given cut¬ 
off value, each angular momentum I from 0 to 12 and parity 
sector (—1)^ and (—1)^^^ contributes as a point in the figure: 
The fact that the points (in red) are superimposed and can¬ 
not be distinguished shows that flmin does not depend on I 
nor on the parity. Furthermore, flmin is always positive, it 
is linear in l/^max a-nd extrapolates to zero for infinite cut¬ 
off (see blue line): This is perfectly consistent with the fact 
that flmin corresponds to the lower border of the 2 + 1 con¬ 
tinuum A.L^i(ik,a), where k has a minimal, discrete value 
scaling as l/xmax in presence of the numerical cut-off, and 
AL=i{ik, ac{2; 1)) vanishes quadratically at fc = 0. In other 
words, there is no negative flmin and no four-body Efimov 
effect for 2 + 2 fermions. 


Appendix A: Is the spectrum of M{s) bounded ? 


In the parity sector (—1)^, the third contribution in 
Eq. (l 6 ^ diverges when {x,u) —>• (0,-1) or {x',u') —>• 
(0,-1). The question is to know if this makes the oper¬ 
ator M{s) unbounded, for a purely imaginary s = iS. 

To investigate this problem, we construct a simplified 
functional that focuses on the diverging part of the ma¬ 
trix kernel ()69l) , replacing each non-zero limit expression 
by its limit, and replacing the vanishing expressions by 
their leading order (here quadratic) approximations: 

M + cha; ~ i(a:^ + y^) with y = tt — 0. (Al) 

Dropping numerical factors and other bounded pieces 
(for example the bit raised to the power s, of modulus 
one), we obtain the mean-D functional 


f dxdu f dx'du' _ •s>'‘ix,u)<s>(x',u') _ 

f dxdu\^{x, u)p 


(D) = 


(A2) 

where the integrals are taken over some convenient neigh¬ 
borhood of {x,u) = (0,-1). It is convenient to use the 
angle 6 rather than u = cosO as integration variable, 
which pulls out a Jacobian sin0 ~ y; we absorb it in the 
integral in the denominator of Eq. (IA2I) with the change 
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of function 

$ = (sin0)^/^$(a:, u) (A3) 

A factor (sin 6* sin ~ {yy'Y^'^ remains in the inte¬ 

grand in the numerator. We restrict the integration over 
{x,y) to the upper half y > 0 of the disk + y'^ < 1. 
Then it is natural to move to polar coordinates: 

i.x,y) = {p cos (j), p sin (j)) (A4) 

so that x'^ + y'^ = p^, x'"^ + and (yy')^^^ = 

{pp'Y^‘^{sm(j>sui(j)'Y^'^. The occurrence of the Jacobians 
p and p' in the elements pdp and p'dp' motivates the 
change of variable in the radial integration: 

A = and X'=p' 2 . (A 5 ) 

Then, considering $ as a function of X and (j>, we obtain 


Since the domain of integration over (A, A') is the square 
[0,1]^, '0 runs over [0,7r/2] and, at fixed 0, r runs over 
[O,i?(0)] with 

i?(0) = min (A12) 

\costp sm ip j 

In the integral over r at fixed 0, we again use the Cauchy- 
Schwarz inequality over the interval r € [0, i?(0)] with 
fir) = |$a(rcos0)| and y(r) = |$o(rsin0)|: 


t-m) 

I (ir|<i)a(rcos0)||$Q(rsin0)| < 
0 



1/2 

■ rRM 

/ dr $a(r cos0)p 


/ dr$a(rsin0)p 

Jo 


Jo 


(A13) 


(fl) = 


fo AW /J''^</>#'(sin0sin0')i/2d)*(A, 0)$(A', 0') 


2/;dA/;d0|$(A,0)|2 

To get rid of the polar angle 0, we introduce 

$„(A)= [ d0(sm0)i/2^(A,0) 
^0 


(A6) 

(A7) 


so that the integral over 0 and 0' in the numerator 
of Eq. (IA6I1 reduces to the product $*(A)$a(A'). In 
that numerator, we use the fact that the modulus of 
the integral over A and A' is less than the integral 

of the modulus, and that ^ ^ ^ ;^/ 2 )i /2 ■ 

In the denominator of Eq. (IA6I) . at fixed A, we ap¬ 
ply over the interval 0 € [0,7r] the Cauchy-Schwarz in¬ 
equality |(/|y)P < (/I/) (y|y) (in Dirac’s notation) with 
/(0) = (sin0)^/^ and y(0) = l>(A, 0); after integration 
of the resulting inequality over A, we get: 


dA|$,(A)|2<2 / dX d0|$(A,0)|2 (A8) 


whose right-hand side is the denominator of Eq. (EH). 
We arrive at 


In the hrst factor in the right-hand side of Eq. (IA13I) , we 
perform the change of variable A = rcos0, so that 

1 pR{' iIj) cos i/j 

/ dr|$a(r cos0)p = - - dA|$Q(A)p 

Jo '"O® r Jo 

where we used i?(0) cos 0 < 1 and the non-negativeness 
of |$aP- The last integral in Eq. (IA14I) is nothing but 
the denominator in the right-hand side of Eq. (lAllIl ! We 
proceed similarly in the second factor in the right-hand 
side of Eq. (IA13I) . except that cos 0 is replaced with sin 0. 
Finally the denominator in Eq. (|Alip cancels out, so that 

PIP) 

1(^)1 ^ / 7 - / / u /2 < +«« ( A 15 ) 

Jq (cos0sm0)^/'^ 

and the spectrum of M{s) is bounded, when s € zK. 


Appendix B: Enunciating and testing a conjecture 
for 64 


I (Si) I < V 


dx |•1■„(V)|2 

We again move to polar coordinates 

(A, A') = (rcos0,rsin0) 


so as to simplify the factor 


1 


(A2-tA'2)i/2 r 


(AlO) 


= — with the 


Jacobian and to obtain 

dr\^a{r cos0)||$a(?'sin0)| 


IMI< 


/o'dA|cI>,(A)|2 


1 . The cluster expansion 


Consider a spatially uniform spin 1/2 Fermi gas at 
thermal equilibrium in the grand canonical ensemble in 
the thermodynamic limit, with a temperature T, and a 
single chemical potential p since the gas is unpolarized. 
The well-known cluster expansion is a series expansion of 
its pressure in powers of the fugacity 2 : = exp(/3y) in the 
non-degenerate limit p —>■ — 00 for a fixed temperature 
T, with /3 = l/iksT) [i^. For our gas, it is generally 
written as 


P\^ 

kYf 


n>l 


(Bl) 


(All) 
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where the overall factor 2 accounts for the number of spin deviations (w) of (w) from the ideal Fermi 

components and A is the thermal de Broglie wavelength gas value [t^ , one has 


/ 2Trh^ A 
yrnksT ) 


(B2) 


When reexpanded in terms of the small degeneracy pa¬ 
rameter pA^, where p is the total density, the cluster 
expansion gives rise to the virial expansion with virial 
coefficients a„ [1^. In practice, one rather considers the 
deviation A 6 „ of from its ideal Fermi gas value, that 
is (for n > 1) from the mere effect of Fermi statistics: 

bn =-^75-1- Abn (B3) 

n°' 

While the cluster expansion has been studied for a long 
time and the second cluster coefficient &2 was obtained 
analytically in reference [i^ (note that bi = 1 accord¬ 
ing to the ideal gas law), there is a renewed interest in 
the cluster coefficients for n > 2. First, the new chal¬ 
lenge is to calculate the for resonant s-wave interac¬ 
tions (with a scattering length a much larger in abso¬ 
lute value than the interaction range), whereas previous 
studies were concentrating on the hard sphere model . 
Second, the bn have been extracted up to n = 4 in the 
unitary limit, from a measurement of the equation of 
state of ultracold atomic Fermi gases [lO, [2l|. The two 
independent groups have reported consistent values of 
the fourth cluster coefficient: 


= 0.096(15) and Abf^'^ = 0.096(10) (B4) 


2. In the unitary limit 

For zero-range interactions with infinite s-wave scat¬ 
tering length a~^ = 0 , i.e. in the unitary limit, the har¬ 
monic regulator method used in (25l |. that introduces an 
isotropic harmonic trapping potential, is quite efficient, 
due to the SO(2,1) dynamical symmetry resulting from 
the scale invariance [H, and the subsequent separa¬ 
bility of Schrodinger’s equation in hyperspherical coordi¬ 
nates in the trap. The value of can be deduced 

from the canonical partition functions, that is from the 
energy spectra, of all the possible fc-body problems in the 
trap, with k < n. One has the following expansion of the 
grand potential of the thermal equilibrium gas in the 
trap: 


-n 

ksTZi 




Br 






(B5) 


where Zi is the canonical partition function for one parti¬ 
cle in the trap, and it is convenient at this stage to be gen¬ 
eral and introduce independent chemical potentials pa for 
the various spin components cr, so that z^- = exp(/3pcr)- 
Then, from the asymptotically exact local density ap¬ 
proximation (see also [H), and introducing also the 


n—1 

2Abn = Y. AS„^.„^=„_„^(0+) (B 6 ) 


where Ai?(0+) = lim^_>Q+ AB{iS) and where we could 
restrict the sum to nc 7 ^ 0 , cr =t, 4 -, since the fully po¬ 
larized configurations are non-interacting and have zero 
deviations from the ideal gas. 

For n = 3, extending to fermions the technique ini¬ 
tially developed for bosons 1 ^. the following analytical 
expression was obtained |29l| |75j| : 


AB2.i(0+)=^ 




— S-^[\nAi{iS,a)] 

7T CiO 


(B7) 

where the function A; is given by Eq. (dl), and the mass 
ratio between the opposite spin component a is equal to 
one (so that Ai? 2 ,i = ABi^ 2 )- It gives 


A 63 ~ -0.355103 (B 8 ) 

in agreement with previous numerical studies and 

with the experimental values [ 2 ^ . 

For n = 4, the problem is still open. A numerical 
attempt [s^, with brute force calculation of the 4-body 
unitary spectrum in the trap, has produced the value 

A^Biume ^ _0.016(4) (B9) 

The disagreement with the experimental results (IB4I1 is 

attributed to uncertainties in extrapolating to w —>■ 0 
the numerical values of (w), in practice obtain¬ 

able only for hui > ksT. An approximate diagrammatic 
theory |24l | (keeping even in the unitary limit only the 
diagrams that have leading contribution in the perturba¬ 
tive regime of a large effective range or a small scattering 
length) gives an estimate closer to the experimental val¬ 
ues (jB4l) . 

Afc^evinsen _ q gg 

Extending the analytical method of reference [ 2 ^ to the 
fermionic four-body problem is technically challenging 
and goes beyond the scope of the present work. On the 
contrary, it is reasonable here to propose and test a guess 
by direct transposition of Eq. (iBTl) : the transcendental 
function A;(s) of the three-body problem is formally re¬ 
placed by det M(^)(s) for the four-body problem, where 
det is the determinant and the operator acting 

on the spinor functions as in the right-hand 

side of Eq. (j 68 p . was introduced and spectrally discussed 
in section IIV BI for the 2-1-2 fermionic problem, and has 
a known equivalent for the 3-1-1 fermionic problem, see 
Eq. (14) of reference [l^. Indeed, in both cases, the 
scaling exponents s (purely imaginary in the efimovian 
channels, real otherwise) allowed by Schrodinger’s equa¬ 
tion in the unitary Wigner-Bethe-Peierls model are such 
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that A;(s) = 0 for n = 3, or such that Eq. (l68t has 
a non-zero solution that is M^^\s) admits a 

zero eigenvalue. Hence our conjecture: 

feN ^ 

r+a2 7 0 7 

X —S—[\ndetMi%^{iS)] (Bll) 

with (n-t.,n 4 ,) € {(1, 3), (2, 2), (3,1)} and is the op¬ 

erator for the four-body problem with ria particles 
in each spin component a. 


3. Existence of the logarithmic derivative of the 
determinant 

The conjecture (IB11|) is not as innocent as it may look 
at first sight. The difficulty is that is actually an op¬ 
erator, and not a finite size matrix: it has a continuous 
spectrum, constituting an infinite, dense set of “eigen¬ 
values” ; even its discrete spectrum may present accumu¬ 
lation points, leading to an infinite but countable num¬ 
ber of eigenvalues. In other words, the determinant of 
is not finite. Numerically, as we have already 
done in section |Vl one of course truncates the unbounded 
variable x to the compact interval [—ccmax, a^max], which 
amounts to imposing the boundary conditions to the 
spinor: 

{x = ±Xi„ax,w) = 0 Vu G [-1, l],Vm^ G {-£, 

(B12) 

After discretization of the x and u variables, M^^^iS) is 
then replaced by a matrix, with a well-defined determi¬ 
nant; still it remains to know if there is convergence of 
the integrand in Eq. (IBllI) when Xmax -boo. As we 
now see, the answer is positive. 

The key point is that what appears in the integrand of 
Eq. (IBllI) is not the determinant itself, rather its loga¬ 
rithmic derivative, which can be written as 

^lndetM(^)(f5') = Tr 

ClD I (ID J 

(B13) 

where Tr is the trace and M ^ the inverse of M. 

Parity (—11^+^: In the parity sector (—1)^+^, the spec¬ 
trum of is at non-zero distance from 0 for a 

mass ratio a = I, as there is no 4-body Efimov effect, 
see Fig. [3l So the inverse of M^^\iS) is well defined. 
Also the operator [iS) is local in the x basis, mean¬ 
ing that the off-diagonal matrix elements of the opera¬ 
tor p-i/ 2 xOTp-i /2 are rapidly decreasing functions of 
|x —x'l, for example there exists a constant such that 

\{x,u,£,Tn^\K^^\iS)\x',u',£,m'^)\ ^ 

[d(x,u)d(x',“ ch(x — x') 

for all x,x\u,u' and all of parity opposite to £, 

and for all .S' G M. Here, we have used Dirac’s notation 


and singled out as in Eq. (1681) a diagonal part and a kernel 
part, 

=V + K^^\iS) (BIS) 

where the operator T) is positive and defined by the 
diagonal-element function d{x,u), 


(D|x, M, £, niz) = d(x, u) |x, u, £, iriz) with 


\ “ 

f 1 H- ^ 1 

e ^ -b ae^ 

_(I + q;)2 

V chx/ 

2(a-bl)chx_ 


(B16) 

This locality is apparent for the first two contributions 
in the right-hand side of Eq. dsg): each contribution is 
bounded, and is consistent with Eq. (jBI4|) at the four 
infinities (x,x') = (±oo,±oo) (see reasoning in section 
IIV Bl) . We then expect that the inverse of that 

can be written as 

= V-^ + (B17) 

is also local, from the geometric series expansion: 

[V + + p-i/2^p-i/2)-ip-i/2 

= + p-l/2 ^('_;L)np-l/2^p-l/2^np-l/2 

n>l 

(B18) 

each term of the series being local (for simplicity, we omit 
to write the exponent (£) and the argument iS). This 
holds of course if the operator is small 

enough. For ^ = I in the 2-1-2 fermionic problem, this 
can be made rigorous: the best constant in Eq. (IB14I) is 

, 2 ( 2 - 

A = — -~ 0.05686 (B19) 

37r 

Then [zl 


l(x,u,£ = l,mz 


= 0lK,(^-^\zS)lx',u',£=l,mi = 0}l 
[d{x, u)d{x', 


< 


2A sh(^|x-x'|) 
y/1 - (27rA)2 sh(2|x-x'|) 


(B20) 


with 6 = arccos(—27rA) G]7r/2,7r[. 

This locality per se is not enough to ensure the conver¬ 
gence of the trace in Eq. (IB13I) . Making the trace explicit 
in that equation and injecting a closure relation leads to 
the writing 


jL 

dS 


Indet M^^^(zS') = f dxdx' f dudu' 

JK. J-1 , 


(x, u, £, rriz \ (i5')] ^ |x', u', £, m(,) 

X {x',u',£,m'J-^M^^\iS)\x,u,£,mz) (B21) 

clD 
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where the sum is restricted to rriz and of parity oppo¬ 
site to that of £ as the exponent (—indicates. The 
locality of and even of -^M{iS), exponentially 

bounds the excursion of |a; — a;'| in the integral over x', 
but still there remains the integral over the unbounded 
variable x. One must take advantage of the structure of 
K(iS) and of its derivative: splitting 

K{iS)= Ki (iS) + K 2 (iS), (B22) 


where Ki and K 2 respectively correspond to the first 
term and the second term in the right-hand side of 
Eq. (l69l) . one finds 

■^K{iS)=t[D,,K,{tS)]+i[D2,K2{iS)] (B23) 


where [A, B] = AB — BA is the commutator of two op¬ 
erators and the diagonal operators Dj are defined by the 
following diagonal functions 


di{x) = 
d2{x) = 




2 cha; 


2 cha; 


(B24) 

(B25) 


such that 


Dj\x, u, rriz) = dj (a:) |x, u, i, m^) (B26) 

Clearly the diagonal term 'D~^ in M~^, see Eq. (IB17I) . 
has a zero contribution to the trace, as [D,Dj] = 0. 
Eq. (jB21() is correspondingly rewritten as 


— Indet = J dxdx' J dudu' 

(x, u, I, m^\K'>2,{iS)\x', u', £, m'^) 

2 

X i[dj{x') — dj{x)]{x', u', £, m'^\Kj^^ (i5')|x, at, m^) 
a=i 

(B27) 

Then di(x) tends exponentially rapidly to 0 when x —>• 
-|-oo, whereas it diverges linearly with x when x —> — 00; 
the contrary holds for d 2 {x). A second property is that 
there exists a constant B such that 


\{x,u,i,mz\Ki{iS)\x',u',£,m'^)\ < 


^+x+x' 


1/4 


B 


dchxchx'/ ch(x — x') 


(B28) 


\{x,u,£,mz\K 2 {iS)\x',u',i,m'^)\ < 

f e— \ B 
Idchxchx'y ch(x — x') 


This is due to the fact, evident from Eqs. (I59I60I) . that 
the denominator in the integral over (j) in Eq. (I69|) . is 
always larger than ch^x — x'). Then, for |x — 

x'\ = 0(1), the upper bound for the matrix elements of 
Ki (respectively K 2 ) tends exponentially fast to 0 when 
X — >■ — 00 (respectively x —>■ -boo), due to the first factor 
in Eqs. (IB28IB29p . which suppresses the linear divergence 
in di{x) (respectively in d2(x)). Then the integral over 
X and x' in the trace converges exponentially rapidly at 
infinity, the logarithmic derivative of the determinant of 
in the (—parity channel is well defined lz3 
and its value can be calculated with a rapidly vanishing 
error in the truncation Xmax when Xmax —>■ +00. The 
numerics agree with this conclusion, and indicate that 
the surprisingly low value Xmax = 5 is sufficient. 

Parity (—1)^: The situation is physically quite differ¬ 
ent for the (—1)^ parity sector, at least for even £: The 
third contribution in Eq. (IMl) is non zero, and it leads 
to a continuous part in the spectrum of that 

reaches zero for even £, see Eq. (15^ . Then the spectrum 
of the inverse is no longer bounded, and its 

matrix elements are not bounded even if one uses the op¬ 
timal (/, tp) representation in which the matrix elements 
of are bounded, see Eq. (IMl), when the lower 

cut-off fmin On the t variable tends to —00. Then, as we 
shall see, there is no exponential locality in the t-basis 
but still the logarithmic derivative of the determinant of 
has a finite limit when tnUn —00, which is 
approached with an error vanishing linearly with 1/tmin- 
To derive this property, a spectral or “Fourier-space” 
analysis is more appropriate than the “real”-space analy¬ 
sis of the previous parity case. After the gauge transform 
and the change of functions performed in Eqs. (I86I89I91L 
the asymptotic t ——00 part of the eigenstates of 
of the continuum (IM)) can be written as 

0^,(t) = fc>0 (B30) 

see Fig. [51 The plane waves and are indeed two 
linearly independent solutions of the eigenvalue problem 
(1921) with 





1 

ch(fc7r/2) 


(B31) 


[see Eq. (15^ . here I is even]. The right solution is some 
specific superposition of these two degenerate solutions, 
with a relative amplitude determined by the physics at 
t = 0(1), that is for (x, u) not extremely close to (0, —1). 
Analytically the value of this relative amplitude is an un¬ 
known function of k and S', but we know that it must be 
of modulus one, so that we can express it as in Eq. (IB30I) 
in terms of a mere phase shift 0(fc, S) € K: (i) the “Hamil¬ 
tonian” for the spinor is hermitian, so that 

the corresponding evolution operator is unitary and con¬ 
serves probability, (ii) the third continuum (IMl) is not de¬ 
generate with the other continua for the considered mass 
ratio a = 1, so the wave incoming from t = —00 
has no channel to escape and must fully get out by the 
incoming channel. 
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gS(t) 

If 

a/ 


exp(ikt) 


a/ 

a/ 


-> 




<- 




-e'*^<^'^>exp(-ikt) 




T- 


A? 

-> 


FIG. 6: Diagram giving the structure of the eigenstates of the 
third continuum (1 84 II . in terms of the variable t of Eq. (EJ: 
fc > 0 is the wave vector of the incoming wave, —k the one of 
the reflected wave with a phase shift 9{k,S). The reflection 
due to the physics at t = 0(1) is, in a toy model, represented 
by a Dirac scattering potential at t = 0 and a hard wall acting 
as a mirror at t = L; the toy model exemplifies the expected 
low k behavior (IB32II of d{k, S). 


The key property that we shall use is that, as it is com¬ 
mon in one-dimensional scattering problems, the phase 
shift 9{k, S) vanishes linearly at low k: 

9{k,S) = kb{S)+o{k) (B32) 

fc—>0 

where b{S) is a S'-dependent effective scattering length. 
We present two plausible arguments to establish this. 
The first argument results from the assumption that the 
writing (jB30|) can be smoothly extended from fc > 0 to 
k < 0: This implies that if one directly replaces k hy —k 
in Eq. (IB30|) . the resulting wave t >->■ 
must reproduce the physical solution ^-ikt 

up to a global phase factor, so that 
and there exists an integer q such that 

-9{-k,S) = 9{k,S)+2qTT (B33) 

The fact that, in one dimension, the arbitrarily low en¬ 
ergy waves are generically fully retroreflected (no matter 
how small, but non-zero, the scattering potential is) leads 
to 9{k, S') —>■ 0 for fc —>• 0 and to g = 0 in Eq. (IB33I) : then 
if 0(fc, S) is a smooth function of k, Eq. (IB32I) holds. The 
second argument utilizes some model for the scattering 
potential in the region t = 0(1), introducing on purpose 
most singular potentials, see Fig. [51 a pointlike fixed scat¬ 
tering center of coupling constant g = (2meffa) placed 
at t = 0 at a distance L from a hard wall (the second el¬ 
ement acts as a mirror and ensures that the wave is fully 
reflected at all energies). At low k, the dispersion re¬ 
lation rifc can be quadratized, 11^ Ri k?k^/{2mes) with 
ruefi > 0, leading to an effective Schrodinger equation 
and scattering problem, so that at fixed S 


Then the property (|B32I) leads to the conclusion that 
the logarithmic derivative of the determinant of 
has a finite limit when the lower cut-off value tmin tends 
to —oo, as we now see. Similarly to Eq. (IB12I) . this lower 
cut-off corresponds to the boundary condition 

= 0 (B35) 

which, considering (IB30I) . leads to the quantization con¬ 
dition for k [T^: 


2fc|tmin| + 9{k, S) = 2n7r, Vn G N* 


(B36) 


Then the contribution of the corresponding eigenvalues 
to the logarithmic derivative of the determinant is 

li, det M<')(iS) I a = y: A In n. = y: ny " 


dS 


n>0 


n>0 


dk dS 

(B37) 

Taking the derivative of Eq. (IB36I) with respect to S at 
fixed n we obtain 


dk 

dS 


ds9{k,S) 


2|tmin| + dk9ik, S) 


(B38) 


In the large |tmin| limit, one can neglect dk9{k, S) in the 
denominator and one can replace in Eq. (IB37I) the sum 
over n by an integral f dn. According to Eq. (IB36I) . 


2—If . I 

^ , rminl ^ 
dn, ^min 


—27r 

— CO 


(B39) 


so that 


— Indet A/(^)(iS')|Q 


r+°°dk 1 dnk 


ds9{k, S) 


dS tmin-^-oo jg 27r rik dk 

(B40) 

This is finite even if flfc vanishes quadratically in fc = 0, 
i.e. it is saved from a naively expected logarithmic diver¬ 
gence, because the phase shift 9{k, S) vanishes linearly 
with k, and so does its derivative with respect to S. 

This main point being established, there remains a 
problem of practical interest, the speed of the conver¬ 
gence with |fmin|- The answer is provided by Poisson’s 
summing formula: 


/(•^”) = T X! 


(B41) 


n^Z 


n^’L 


for any A > 0 and for an arbitrary function /(fc), f{x) = 
Jg^dkexp{—ikx)f{k) being its Fourier transform. For 
simplicity, we give details in the case where 9{k, S) is 
linear in fc at all fc, that is 9{k,S) = kb{S). From the 
quantization condition (IB36I) one has 


^ (fcg) ^-b(tanfcL) ^+i 

{ka)~^ + {tankL)~^ — i 


(B34) 


The phase shift 9{k, S) is indeed an odd function of fc, 
and at low fc one indeed obtains the linear law (|B32I) with 
2/6(5) = l/a{S) + l/L{S). 


so that 


27r 

k = Xn with A = —;-;-- 7-^7 

2|fmin|+6(5) 


dk db{S) A 
dS ~ d^^^ 


(B42) 


(B43) 
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This, together with Eqs. (jB37l) . leads to a function / 
given by 


f{k) 


J, db{S) 

^ dS 


drik 

dk 


such that 

^ Indet ^ /(An) 

n>0 


(B44) 


(B45) 


Then using the fact that the function / is even, one can 
express the sum over N* in terms of the sum over Z, and 
then in terms of /(O) and /: 


n>0 nGZ 

= -^/(0) + |^EA2WA) (B46) 

n£Z 

The function / is a smooth function of fc, in particular 
in /c = 0, that rapidly decreases at infinity, so that its 
Fourier transform f{q) is also rapidly decreasing when 
l^l —?► +00. In the large /mini limit, 1/A diverges linearly 
in /mini and one commits an exponentially small error 
0[exp(—C/minl)] (C is some constant) in neglecting the 
n 0 terms in the last sum over n in Eq. (IB46I) . As a 
consequence, 


jL 

dS 


1 

Indet Mi^i(iS')|0 = dk 

^min ^ OCi STT Jq 


+ 


db{S) 

dS 


2/n 


b{S) 


fife dk 
+ 0[exp(-C'/min|)] (B47) 


4. Other convergence issues 


To show that the conjectured values (IBllI) are finite, 
one must also check that the integral over S is conver¬ 
gent at infinity, and that the sum over the angular mo¬ 
menta £ is convergent. This we have initially explored 
numerically. First, at a given i, we found that the loga¬ 
rithmic derivative of the determinant of {iS) rapidly 
decreases when S —>■ -l-oo, presumably exponentially fast, 
see Fig. [ZKa). Second, after the integration over S is 
taken, one observes also a rapid convergence of the series 
over £, see Fig. EKb), if one takes the precaution to be 
accurate enough in the discretization of the integral over 
u Is^l. 


These numerical results suggest that the contribu¬ 
tion of angular momentum i to Ai32™'^(0+) and to 
Ai?3°/'’(0+) can be obtained, when £ is large enough, 

from a perturbative calculation in Eq. (IB13I) . lim- 

ii) 

ited to leading order in the operators A'j ' defined by 
Eqs. (IB15IB22I) . at least for the 2-1-2 problem in the par¬ 
ity channel (—1)^“'"^ where ATinv in Eq. (IB 171) has a chance 
of being bounded. This idea was implemented with suc¬ 
cess at the three-body level in reference , treating the 
integral term in Eq. (l78l) as a perturbation of the constant 
term [8l|. 


Let us implement the idea for the 2-1-2 problem, in the 
parity sector (—1)^+^ of the subspace of angular momen¬ 
tum £. We truncate Eq. (IB18|) to order one included in 
the operator K, to obtain 


d 

— In det M ~ Tr 
db 


{p-^ - V-^KV-^) -^K 

CltJ 


(B49) 


where we have replaced A, /(O) and /(O) by their values. 
When 9{k, S) is not a linear function of k, we obtain the 
general result 


d 1 

— lndetM(^i(iS')|0 = dk 

C/O ^min ^ OC) 2i7T Jq 

kd^dsOikp) 


ds9{k, S) dOfe 


-I- lim 


fc—>0 flfe[2/n 


dkOikp)] 


fife 

0[exp(-C'/n 


dk 


(B48) 


In any case, when —oo, the limiting value of the 

logarithmic derivative of the determinant of is 

approached with an error that vanishes only polynomially 
with l/tmin [Zl 

On the contrary, if the dispersion relation fife nowhere 
approaches zero, as for an odd £ in the parity sec¬ 
tor (—1)^, the limfe_).o term in the right-hand side of 
Eq. (IB48I) is zero and the convergence of the logarith¬ 
mic derivative of the determinant is exponentially fast 
with /mini, as also observed numerically; this last situa¬ 
tion is then similar to the exponentially fast convergence 
of ^lndetM(^)(i5') when Xmax +oo, which is always 
achieved for the continua (1771511) . even when the third 
term in Eq. (15^ is active. 


Then we split K as in Eq. (IB22() and we use the commu¬ 
tator structure (IB23I) . Using the invariance of the trace 
in a cyclic permutation and the fact that the diagonal 
operators T> of Eq. (IB16I) and Dj of Eq. (IB26I) commute, 
we find that only the crossed quadratic contributions in 
Ki and K 2 survive, so that 


dS 


In det M ~ Tr 


-p-i/TiH-i—- (1 O 2) 


= ^ Tr {-V-^KiV-^K^) (B50) 


Integrating by parts in Eq. (|Bllll and using the fact that 
the integrand is an even function of S we obtain the ap¬ 
proximation 


Jq tt dS Jr2tt ^ 

(B5I) 

Calculating the trace in the \x,u,£,mz) basis (with I -|- 
TTiz odd) and injecting a closure relation in that basis as 
e.g. in Eq. (IB21|) . we realize that the integrand has a 
very simple dependence with S', due to simplifications as 
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FIG. 7: (Color online) Convergence in the integration over S and summation over i in Eq. (IBllI) (we recall that the mass 
ratio is a = 1). (a) The logarithmic derivative of the determinant of is a rapidly decreasing function of S', the figure 

takes as an example (al) the £ = 0 channel of the 2 + 2 fermionic problem (for a numerical cut-off fmin = —9, thus without 
extrapolation; note the minus sign in the vertical axis), (a2) the i = 0 channel of the 3 + 1 fermionic problem, (b) The sum 
over £ also seems to converge well, see the contribution of each angular momentum channel to the result (IBllll for (bl) the 
2 + 2 problem and (b2) the 3 + 1 problem: black disks for the parity sector (—1)^, red disks for the parity sector ( — 1)^^^; the 
plus signs indicate the corresponding cumulative sums. Convincing evidence is even shown in (b3) for the 2 + 2 problem in 
the (—1)^"*"^ parity sector and in (b4) for the 3 + 1 problem in both parity sectors, where the numerical results (black disks for 
parity (—1)^, red disks for parity (—1)^^^) are compared to the perturbative results (IB54II and (IB60II [red asterisks for parity 
( —1)^, black asterisks for parity (—1)^^^] that extend to four bodies a technique developed for three bodies in reference 
and are expected to be exact asymptotic equivalents for £ —>■ +oo (what is actually plotted is the absolute value of the results, 
so as to allow for a log scale, but their sign is indicated with the label “< 0” of the same color as the corresponding disks 
when they are negative: the negative black (red) disks are indicated with a black (red) “< 0” label; note that the black (red) 
asterisks have always the same sign as the corresponding red (black) disks). 


follows: 


/ ch x' \ 
chx ) 


s/2 



cha; 

cha;' 


s/2 


(B52) 


decreasing function of x or x', so that, as a function of 
S, it is a rapidly decreasing function. This gives a reason 
for the numerically observed fast decay of ^ In det M at 
large S. Also, integration over S is straightforward due 
to 


where we wrote the phase factor of the first term as it 
is in Eq. ((69ll and the phase factor of the second term 
of Eq. (l69ll with x x', and used s = iS with S real. 
So the integral over x or x' takes the form of a Fourier 
transform with respect to x or a;', with S as the conjugate 
variable; this is the Fourier transform of a smooth rapidly 


f ^i(x — x')S 
JR 27r 


S(x — x) 


(B53) 


We finally obtain the leading order approximation 
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2,2 parity( 1)^+1 Jq Jq Jq 2tT J^ 27 ^ d{x,u)d{x,u') c\\X 

J E E (0 + nTT, O' + n'TT, 0') 


n—Q n'—{) 


{ 1 + [(u + e-^){u' + e-") + W cos<^l 1 11 + [{u + e^){u' + e") + vv' cos(j>'] \ 


(B54) 


In the integrand of Eq. (IB54E d{x,u) is given by Eq. (IB16E we again use the notations u = cos 9 and v = sin0 and 
the same for 0', and we introduced the function 




mz,m'——t 


= Tr, 


^—i9Lz/h^i<j>Lx/h^i9'Lz/h^i<j>'Lx/h 


sin[(2^+ 1)^2] 
sin((5/2) 


(B55) 

(B56) 


where the trace is taken over the whole subspa ce U l^rriy), —I < mz < of angular momentum I without any parity 
restriction and the angle 5 £ [0,7r] is such that 


1 + 2 cos d 


uu{l + cos (j) cos (j)') — {u + u) sin^sint/)' + ur;'(cos </> + cos cj)') + cos (j) cos (j)' 


(B57) 


The sum over n and n' in the numerator of the inte¬ 
grand of Eq. (IB54I) suppresses the contribution to Tt 
of the states \^^mz) and |^,m^) of the wrong parity, 
(_l)m^ _ (—= (—1)^. We expect the approxi¬ 
mation (|B54|) to be an exact asymptotic equivalent for 
^ —>• +00, and this is also what the comparison to the 
numerical results in Eig. [7l^b3) indicates. Amazingly it is 
already good for ^ = 1, as it deviates from the numerical 
value by about 9% only. 

This perturbative treatment can also be applied to the 
3 + 1 problem, using the integral equations of reference 
The main difference is that the spinor 
is now subjected to a condition reflecting the fermionic 
exchange symmetry of the two j, particles that are spec¬ 
tators of the interacting pair [Ig , 

4>E. y) = (-l)^+i$W (a,, u) (B58) 


This means that the kernel must be restricted to 
the corresponding subspace, hence the occurrence of a 
projector P = (1 + [/)/2 on that subspace, where the 
unitary operator U such that in Dirac’s notation 

U\x,u,i,mz) = - x,u,£,-mz) 

= -P^^‘’’/^\-x,u,£,m^) (B59) 

is an involution (t/^ = 1) [s^. The interesting point is 

now that, even if -^K is a sum of commutators as in 
_ _ I dS 

Eq. (IB23p . the corresponding Dj do not commute with 
the projector P. As a consequence, when one expands 
M~^ up to first order in K, ^IndetM contains both 
a contribution of order one in K and two contributions 
of order two in K. Here is the result in the subspace of 
angular momentum £ and parity e [s^: 
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1 1 

- ^ e”7i(0 + rnr, 0, </> + tt, 0) 


27tV2Jo Jo 2^ rf3i(0,«) 3^ J^(2 u + m2 + ^;2cos</.) 

I + a 


^-^fd^Tdordo' 


vv 


JO ^0 Jo 27r Jo 27 r ^3,1(3;, u)d3,i(a;, m') dia; 
1 ^ ^ 

- ^ ^ e”+” T«(6» + rnr, 9' + nV, (jj, 4>') 

n—0 n'—O 


|2 + e [e ^{u + u') + mm' + mm'cosi^] | 12 + [^^{u + m') + mm' + mm' cos (/)']| 


^-^jd.' rdo rdo' 


j \ 1/2 


'0 Jo 
1 ^ ' 


27r Jo 27r Icha;' / ^31(0,M)(i3i(a;', m') 


- ^ ^ e:”+” (0 + MTT, m'tt, </>, </)' + tt) 


n=0 n'—O 


2e“ 


2a 


+ + --(ue “ + m' + mm' + mm' cos 

1 + a 


, , 2a 

2e““ + + :;-(Me““ + m'+ mm' + mm' cos </>') 

1 + a 


(B60) 


where d^\(x, u) defines the diagonal part V of the opera¬ 
tor M for the 3-1-1 problem (as d{x, u) did for the 2-1-2 
problem), see reference [T^ : 


d3i(a;,M) 


1 -f 2a au 

(1 -I- a)^ (1 -I- a)2 chcc 


(B61) 


As can be checked in Fig. [T^bd), this approximation is 
in good agreement with the numerical results even for 
i = 0, where it deviates from the exact result only by 
~ 13%. In the large i limit the first contribution in the 
right-hand side of Eq. (IBBOp rapidly dominates over the 
other two; summing over the two parity sectors e = ±1 
and restricting for simplicity to a mass ratio a = 1, one 
can integrate it over 9 and (j) at fixed S € [0,7r], where 
1-1-2 cos S = u + cos 4> + u cos cj) as shown by Eq. (IB57I) 
taken with 0' = </>' = 0, to obtain the rapidly decreasing 
large-£ equivalent (s^ 




2 i+l 

^->-1-00 27r2 


dS sin[(£ -I- 1/2)(5] 


arccos■ 


8 cos^ i5 -I- 5 cos i5 — 1 
3(3 -I- cosi5) 


[(5 -I- 4cos(5)(1 -I- cosh -I- cos^ h)]^/^ 


(B62) 



FIG. 8: (Color online) Fourth order cluster coefficient ()B67I) 
for a harmonically trapped unpolarized spin 1/2 unitary 
Fermi gas at temperature T, as a function of Phw, with 
P = l/(fcsT), u) the angular oscillation frequency of a fermion 
in the trap and a = mrlmi = 1. Blue line with symbols: re¬ 
sults of reference [s^ obtained by brute force numerical cal¬ 
culation of the up-to-four-body spectra in the trap (disks: 
actually calculated values; circles: values resulting from an 
extrapolation). Red lines: our conjecture (IB 66 II (the values 
slightly differ depending on the linear or cubic extrapolation 
to the numerical cut-off limit 1/tmin —>-0“). 


leading, after use of Eq. (IB6L to 

5. The verdict 


The numerical results for our conjecture (IBlip are 

AB2°2"j(0+) = -0.0617(2) (B63) 

AB3™j(0+) = -b0.02297(4) (B64) 


A64™J = -0.063(1) (B65) 

This clearly disagrees with the experimental values (IB4I) . 
Remarkably, for our conjectured value is very 

close to the approximate diagrammatic result 0.025 of 
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[^ . whereas for AB 2 . 2 { 0 '^) it widely differs from the 
(still approximate) result —0.036 of [2^ (these values 
were communicated to us privately by Jesper Levinsen). 

A useful complementary test is to compare to the theo¬ 
retical results of reference . As mentioned above and 
in that reference, these results, obtained with the har¬ 
monic regulator technique, are trustable at non-zero val¬ 
ues of phu! without extrapolation to = 0 (w is the an¬ 

gular oscillation frequency in the trap and /3 = l/(kBT)). 
It is actually straightforward to extend with the same no¬ 
tations the conjectur e (IBllI) to a non-zero value of ui, see 
Eq. (38) of reference |28l| : 


^eN 

dS smjSptiuj) d 

^ Jq tt sh^phuj) dS ^ 


[lndetMW„^(i,5)] (B66) 


Since | sin(5'/?/ia;)/sh(/3/ia;)| < S, this does not raise 
new convergence issues and the numerical evaluation of 
Ai? 2 ° 2 ”'^(w) and Ai? 3 °"-^(a;) is straightforward once the 
logarithmic derivatives of the determinant of M are 
known. The resulting value of the fourth in-trap clus¬ 
ter coefficient 

ABi{ui) = -[Ai?3 i(w) -b Ai32,2(w) -b Ai?i_3(w)] (B67) 


(with Ai? 3 _i = ABi _3 for the mass ratio a = 1) is plotted 
as a function of PHcu in Fig . [5] It clearly disagrees with 
the results of reference [^, not only with the ones re¬ 
sulting from the extrapolation to phio = 0 but also with 
the actually calculated ones. 

The conjecture is thus invalidated, and more theoret¬ 
ical work is needed to derive the correct analytical ex¬ 
pression for AIi 4 of the unitary Fermi gas. 
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dimensional J(u„ • k) where (u„) is an orthonormal ba¬ 
sis. As explained in the text one can take k = k 2 -|-k 4 -|- 
77 v(/ 3 ) 77 z( 7 )(kf’' + kf’'). We first take ui = ey so that 
ui-k = ( 77 z(- 7 )ev)-(kf-fkf) = sin 7 ex-(kf-fkf) = 
— sin 7 sin/3oIkf-|- kf’'] where we used Eq. (I5IJ1) . This 
gives the factor ^(sinq) in Eq. (1491) . As explained in 
the text we can restrict to the case 7 = 6 (up to a 
change {3 —/3) and we are left with a two-dimensional 

Dirac (5(k_L) in the plane orthogonal to ey. In princi¬ 
ple ky = k 2 -I- k 4 -f 77y(/9)(kf^ -I- kf’') but we can 
equivalently take ky = 77y(—/ 3 )(k 2 -I- k 4 ) -|- kf^ -|- kf’' 
due to the rotational invariance of the Dirac distribu¬ 
tion. Using 77y(—/3)ez = cos fdez — sin/Sex and tak- 

ing U2 = I, L 11 Li = cos/3oez — sin/3oex and its dual 

1^1 H-kg I 

ua = sin /3o ey -|- cos Pc ex in the ZX plane, we justify 
Eq. (l49ll . 

[54] Similarly to Eq. (l53l) . ]£,m^, = 0) = 

\t,mz = 0 ) = mz = 0 ) where sy are 

just signs since \l, = 0 ) can be taken with real compo¬ 

nents in the \i, m^) basis {Ly has purely imaginary ma¬ 
trix elements). Then S-S+{£,mz = m^ = 

0) = 1. The action of jjj Cartesian coordinates 

is {x,y,z) {—x,y,—z); in spherical coordinates of po¬ 

lar axis z it is {0,(f>) —^ {tt — 0,% — p). For YP‘^’^{0,<j)), 
which does not depend on p, this is equivalent to the 
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action of parity { 6 , (p) ^ {n — 9, tt + (p) and it pulls out 
a factor (— 1 )* so that s- = (—l)^s+ and \£,mx = 0 ) = 

= 0). Se¬ 
ries expanding the exponentials in this last expression, 
and using the fact that Ly only couples states of differ¬ 
ent niz parity, one gets Eq. CT . 

[55] One can keep [0, 27r] as the range of integration over (p' 
since the integrand is a periodic function of (p' of period 
27r. 

[56] For 6 o > 6 i > 0, -5i)(5o + 

6 i)]^/^ with zo = -bi/{bo + [( 6 o - bi)(bo + ^i)]^^^}- 

[57] If a; —^ —oo, x' —> -foo or a; —>■ -foo, x' —>■ —oo, the matrix 
kernel (1691) entirely tends exponentially to zero, which 
does not bring any significant new information. 

[58] This is trivial for ^ = 0. For £ > 1, this results from 

the fact that, for L > 1, one can take rrix = 1. Then 
{£,mz\£,mx = 1 ) 7 ^ 0 , except if £ is even and = 

0 (in agreement with Eq. (I54E considering the x y 

symmetry), in which case one may return to the choice 
rrix = 0 and use the fact that (£, = 0 ]£, rrix = 0 ) yf 0 

for even £. 

[59] Let us explain more physically why the A_l function ap¬ 

pears in the expression of the continuous spectrum. The 
idea is to consider a physical state of the 2 - 1-2 fermionic 
system, corresponding to a non-zero square integrable 
solution of Eq. (I 68 E and to see how the four- 

body wavefunction scales when three particles, say 1 , 2 
and 3, converge to the same location, the fourth parti¬ 
cle being at some other fixed location. As we have seen, 
an extended eigenstate of the continuum that varies for 
X ^ -00 as = 0 )Pl{u), k € R, 

has an eigenvalue Q = AL{ik,a). According to the ana¬ 
lytic continuation argument of reference [l^] , this implies 
that the = 0 localized eigenstate "I?ml [x, u) vanishes for 
X —>■ —00 as mz\£, rtix = 0 )Pl{u), where the 

real quantity k is the minimal positive root of 

Al{k,, a) = 0 

with L chosen so as to minimize k (minimizing k amounts 
to selecting the most slowly decreasing exponential func¬ 
tion e'^^, that is the leading contribution for x —^ — 00 ). 
This implies that k is one of the possible scaling expo¬ 
nents S 3 of the 2-1-1 fermionic problem, see Eq. (1821) . In 
order to determine the limit of Eq. (unj when ri 3 and 
]r 2 — R 13 ] —>■ 0 tend to zero with the same scaling law, 
that is both are proportional to the vanishing hyperradius 
Ri 23 of particles 1, 2 and 3, we determine the large k 2 
limit of the integrand of ip 2 A in Eq. (unj at fixed k 4 : Omit¬ 
ting to write the angular part for simplicity, we find that 
D(k 2 ,k 4 ) scales as so that its Fourier 

transform, according to the usual power-law counting ar¬ 
gument, scales as ]r 2 — Ri 3 ]'‘“^]r 4 — Ri 3 ]®“'‘“®^^. The 
same reasoning applies to ipiA. At fixed ]r 4 — R 13 ] > 0, 
the four-body wavefunction therefore scales as Rl^-i = 
exactly as predicted by Eq. (5.179) of reference 
[ 34 ]. This whole discussion is formal for the 2+2 fermionic 
problem since, as we shall see, there is no four-body Efi¬ 
mov effect, but it explicitly applies to the 3 + 1 fermionic 
problem and nicely completes reference [l^ . 

[60] The continuous spectrum can be recovered by 

keeping only but exactly the last contribution in Eq. (1691) 
to the matrix kernel of Eq. (I 68 II , that is without 


resorting to a local approximation of this contribu¬ 
tion around {x,u) — (0,-1). The explicit calculation 
remains simple for a unit mass ratio a = 1. The 
eigenvectors of the resulting operator are then of the 
form ^l^^{x,u) = {£,mz\£,mx = 0)^{x,u) 

with the ansatz >I>(a;,'u) = {chx)~^’‘~^^'^P'^{u + 

cha;)“*'®'''^/^^'^^x('\/2fc24/A'24), ^24 and K 2 A being the 
relative and center-of-mass wave numbers of particles 2 
and 4, so that 2 k 2 AlK 2 A = i chx+u ')^^'^■ then obtains 
the integral equation for x(fc): £l{ 2 k^ + l)^^^x(fc) = 
(fc2 + l)i/2^(fc) _ J+“ dk'k'\ik')/il + e + fc'=). 

Further setting x(^) = k~^^^{l + 2k'^)~^^'^%p{t) with 

/ — 2t \ 

k = exp(—t), one obtains Q.%p{t) = i>{£) ~ 

2(-l)^ r Jj/_ </;({') exp[-3(t-n')/21 _ rp, 

TT Jr“'' (l + e-2t+e-2t')[(l + 2e-2t)(l + 2e-2t')]l/4' 

t ^ —OO continuum of that eigenvalue problem solves 
(%/2n - l)V'oo (t) = /r dt' plane wave 

solutions ipooit') = reproducing (18411 . The ansatz 
for ${x,u) results from the fact that D(k 2 ,k 4 ) is of the 
form ■ K24/A:24)x(V2fc24/A'24) {Pt is a 

Legendre polynomial), which is apparent if one turns 
back to Eq. m and realizes that its last contribution 
conserves the total wave vector K 24 (up to a sign). 


[ 61 ] In reality, the physical solution for the tensor >I>ml(a:,u) 
must correspond to a zero-eigenvalue of the matrix 

In this respect, fc = 0 is acceptable only for 
£ even. Furthermore, as we shall see in Appendix 
[see the footnote called above Eq. (IB36II ]. when £2 = 0, 
(p{t) actually scales linearly in t for t —>■ —00, so that 
D(k 2 ,k 4 ) actually diverges as ln(|k 2 + k 4 |/fc 24 )/(|k 2 + 
k4|/fc24)“'''^'^^ when k 2 + k 4 — >■ 0. This would result 
in a In(r 24 /|R 24 — Ri 3 |)/r 24 divergence of the function 
^I 3 (r 2 — Ri 3 ,r 4 — R 13 ) at r 24 = 0, which is physically 
opaque. Let us keep in mind however that, according to 
section |V] there is no four-body Efimov effect in the 2 + 2 
problem, so that Eq. (1681) does not actually support any 
non-identically zero solution 4?^^^ if s G iR. 

[ 62 ] For £ > 0 , one can use the identity: (e ■ 

62,64x2 • 62,624 • 62)]*e*"‘*^^='‘+®^*/^'(£,m2]£,m2: = 0) = 
s+Y£°(K 24) where K24 = (k2 + k4)/]k2 + k4], we recall 
that 724 = r24 — 024 /2 and r24 is the angle between k2 
and k2 + k4, the spherical harmonics notation is as in 
footnote [ 51 ] and the sign s+ is defined in footnote [ 54 ]. 
This identity is also useful for footnote [ 60 ]. 

[ 63 ] If k24 —>■ +00, ]a:] necessarily tends to +00 since ]u] < 1 , 
in which case K2A and k24 both diverge as and have 
a non-vanishing ratio. 

[ 64 ] We used dysin(ay)sin(by)/y'’~^^''^ = |r(i — s) 
sin[-|(l — 2 s)][]a — fo]'’“2 — ]a + 6]'’“2], for s £ iR, a and 
b real numbers that differ in absolute value. 

[ 65 ] The integrand in Eq. (IWl) is 0 (l/q^) when q —>• +00, so 
that the integral converges when r24 —>■ 0 at fixed non¬ 
zero R24 — R13. Furthermore it converges to a non-zero 
value (e.g. to tt for s = 0). 

[66] If one takes a function D(k2,k4) with no singularity in 
k2 + k4 = 0 , for example ‘I>Q°^(a:,u) = exp(— k]®]) in the 
ansatz (HSJ as in footnote [ 59 ] for £ = L = 0 , one finds 
by an explicit calculation that .4(r2 — Ri3,r4 — R13) is 


finite for r 2 = r 4 = R 24 . 


[67] There is a paradox here. As shown by Eq. @, the ac- 
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tion of the Hamiltonian H on i/) 24 (ri, r 2 , ra, r 4 ) leads to 
a (5(ri — rs) distribution, not to a 5{r2 — r4) distribution. 
This shows that ip 24 is the so-called 1 — 3 Faddeev compo¬ 
nent, and it cannot have any l/r 24 singularity. How can 
yli 3 then have such a singularity? The answer as usual 
lies in the order of the limits. At fixed non-zero ria, it 
is apparent that the function u(ri 3 ) in Eq. (IIUI) . through 
its dependence on gia given by Eq. (HU, provides an ul¬ 
traviolet cut-off of order l/ria in the (k 2 , k 4 ) wavevector 
space, so that tl >24 cannot diverge when r 24 —>■ 0. But if 
one first takes the ri 3 —>■ 0 limit, the function w(ri 3 ) is 
replaced by its equivalent l/(47rri3) which has no mo¬ 
mentum dependence: The wavevector cut-off is set to 
infinity and a l/r 24 divergence in limri 3 ->o(ri 3 '(/> 24 ) can 
now take place at r 24 = 0. 

[68] This reasoning can be transposed to the case of four iden¬ 
tical bosons, when three of them converge to the same 
location in the relative three-body channel where the Efi¬ 
mov effect takes place. As this channel has a zero angular 
momentum and an even parity, this implies that, in such 
a configuration, the total internal angular momentum £ 
of the four-body system is carried by the relative motion 
of the fourth boson with respect to the center of mass of 
the first three bosons, leading to a global parity (—1)^. 
This indicates that, for i ^ 0, the four-boson unitary sys¬ 
tem in an isotropic harmonic trap should have interact¬ 
ing states in the (—1)^"*"^ parity sector that are immune 
to the three-body Efimov effect. Such “universal” states 
have indeed been observed numerically in reference [d^j . 
but for a total internal angular momentum i = O: this 
observation cannot be explained by our reasoning. 

[69] If for simplicity one omits to write the other variables, 
the integral over 9', of the form d9' 

is approximated by '^^^iW{9j)K{9i,9j)${9j), where 
{di)i<i<ng is the set of (non-equispaced) discrete values 
of 9 proposed by the ng-points Gauss-Legendre method 
and w{9i) the corresponding weights. To render the re¬ 
sulting discretised form of M(s = 0) hermitian, it suffices 
to take as unknowns w{9j)^^^^{9j) and to multiply the 
eigenvalue equation by w{9i)^^^, which leads to the ker¬ 
nel [w(9i)w(9j)]^^^K{9i,9j) without modifying the spec¬ 
trum. 

[70] With the notations of the previous footnote, one takes 

[dx (x,9i) as unknowns in the zone p > po, 

and [dt {x{t,il)i),9{t,ij)i)) as unknowns 

in the zone p < po- In this way, after multiplication of the 
eigenvalue equation by \dxwx{di)\^^^ or \dtw{ipi)\^^^ poe^ , 
one obtains an hermitian matrix. The ^-dependence of 
the weight Wx{di) results from the ^-dependence of Pmax. 

[71] The triangle corresponds to the zone H < 

Al=i(0, a). 

[72] The same conclusion must hold for the eigenvalues in 
Fig. [3] [parity (—1)^^^] that are above the upper external 
border of the continuum, as well as for those below the 
lower external border of the continuum (there are some 
for a close to unity). 

[73] It is explicitly supposed here that a potential purely 
imaginary root S 4 of Eq. (I22II would exist above some 
threshold value ac(2; 2) of the mass ratio a, with S 4 = 0 
at threshold. One can however imagine another scenario, 
with S 4 still a continuous function of a: S 4 would exist for 
all mass ratio a G [1, q;c( 2; 1)], with a;c(2,1) = 13.6069 ... 
the three-body Efimov effect threshold, in which case S 4 


would not need to cross zero for some a. This scenario 
is however excluded (i) by the experimental results for 
the spin 1/2 unitary Fermi gas, which has a mass ratio 
a = 1 (no significant four-body losses are observed) and 
(ii) by the numerical calculations in the conjecture on the 
fourth cluster coefficient 64 of that unitary Fermi gas in 
Appendix [B] fit is found for a = 1 that the operator M(s) 
is positive for all purely imaginary s, which excludes the 
existence of a root S 4 ]. 

[74] Note that is actually equal to as soon as 

the two indices differ from zero, since the ideal gas grand 
potential is the sum of the grand potential of each spin 
component, and no t- 4 - crossed term can appear in the 
resulting cluster expansion. 

[75] It can be deduced from Eq. (7) of reference [ 2 ^ by inte¬ 
gration by parts. 

[76] For £ = 1 within the even parity sector, the minimal value 

of A given by Eq. (IB19II corresponds to x = x' ^ -foo, 
u = m' = 0 in Eq. (IB14I) . Let us use the notation 
O ^ {£ = l,mx = = l,mi = 0) 

and introduce the operator T in the space of functions 
of the single variable x such that, in Dirac’s notation, 
(a;|T|a;') = A/ch(a; — x'). Then Eq. (IB14II can be 
rewritten as |(a;,u|0|a;',u')| < (a;|T|a:'). For any integer 
n greater than 1, we inject n — 1 closure relations and 
use the triangular inequality to obtain \{x,u\0'^\x',u')\ < 
J^dxi ... dxn-i dui ... dun-i(x\T\xi)... (a:„-i|r|a;') 
= 2"-i(a:|r"|a:'). Then |(a:,u|[(l -f O)-^ - t]\x',u')\ < 

where we have 

used a series expansion in powers of O and where 
tk = ~ 7 rA/ch(A: 7 r/ 2 ) is the eigenspectrum 

of T. Calculating the integral over k and combining 
the result with the identity {£ = l,mx = [^ = 

l,m(, = 0 ) = r»-i/=[(]l -f O)-^ - ]l]r»-i/= leads to 
Eq. (IB 20 I 1 . 

[77] To finish the calculation, one can use the mild hy¬ 
pothesis that the matrix elements of Ainv are uni- 

/ x + x' \l/4 

formly bounded, that |di(a:) - di( 2 : )| < 

(aThr) (afhT^) , due to exp(a;') < 

2 cha:' or exp(a;) < 2 cha; and to the triangular inequal¬ 
ity |di(a:) — di(a;')| < |di(a;)| -|- |cii(a;')|. Then one can 
integrate over x' or over x (depending on the term), us¬ 
ing Jjj ch(t- 3 :') ~ finally faces the integral 

/ x \i/4 

<+ 00 . 

[78] The value k = 0, that is n = 0, should not be included. If 
one directly takes the limit fc —>■ 0 in Eq. (IB30II one gets 
the absurd result 4ik=o{t) = 0. The correct way of taking 
the limit is to first divide Eq. (IB30II by ik. One then finds 
that (j)k=o{t) diverges as 2t — b{S) when t —>■ — 00 , so it 
does not satisfy the boundary condition (IB35I) . 

[79] In the numerics, we extrapolate to 1/tmin = 0 using 

a cubic fit in 1 /tmin, with data down to minimal val¬ 
ues l/|tmin| = 1/200 for .^ = 0 and l/|fmin| = 1/30 for 
£ > 0. For ^ = 0, as a test of the finite-tmin formalism, 
we have used Eq. (IB48I) to predict the leading numerical 
error on °^( 0 ^) due to the tmin truncation, that 

is ( 87 r|tmin|)“^/g ^^[^(oo) - b(S)] = 2.3(l)/(87r|tmm|), 
which agrees with the direct numerical calculation. To 
obtain b{S) at any given S, and hence the integral of 























30 


6(oo) — b{S), we calculated numerically the eigenvectors 
corresponding to the first few eigenvalues Qn {n > 1) 
of = iS), and we fitted the corresponding func¬ 

tions (jinit) [defined as in Eq. (I91ll 1 with a three-parameter 
sine function t sin(fc„t — 9n/2) as suggested by 

Eq. (IB30I) . where An is a complex amplitude, kn an ef¬ 
fective wavenumber and 6^ a phase shift. The fits are very 
good, the obtained values of kn agree very well with the 
dispersion relation (1841) : setting 6{kn,S) = 9n, we also 
find that the quantization condition (IB36II is well obeyed; 
finally, extrapolating 9n/kn to n = 0 linearly in k^ gives 
b{S). To be complete, we note that b{S) looks like a nega¬ 
tive amplitude Gaussian on a non-zero background fe(oo), 
that is fo(oo) — b{S) ~ 1.05 x exp(—0.6685'^). As the vari¬ 
able t in Eq. m depends on po, so does b{oo). In our 
numerics, po = 2/5 and we find b{oo) = 3.84(1). More 
analytically, one expects at large S that the phase shift 
9{k, S) is imposed by the third contribution in Eq. (I69II . 
the first two ones becoming rapidly oscillating and neg¬ 
ligible. Then one can use the analytical results of foot¬ 
note [60]: introducing the phase shift 6(p(fc) such that 
tp{t) — sin[fct — 9^{k)/2] + o(l) for t —>■ —oo, one expects 
that 9{k,S) —>■ 9^{k) — 2fcln(po/'\/2) for S —>■ oo, so that 
6(oo) — b^ — 21n(po/'\/2) with 6^ ~ 1.33. Our numerics 
fulfill these expectations, which constitutes a good test. 

[80] In practice, we used a Gauss-Legendre scheme with up 

to 59 points, using 9 rather than u = cos 9 as the inte¬ 
gration variable, with the change of function (lA3l) and 
the inclusion of the extra Jacobian (sin 6 sin in the 

matrix kernel. 

[81] The integration over S' G R of Kf,{iS )/ cos u — 1 using the 
first line of Eq. (|78l) leads to duPe(u)/(l -(- usini/). 
The large i limit of that integral reproduces exactly 
Eq. (42) of reference [2^, as we have checked using 
Pe{u) = {2^£\)~^— 1)^], then integrating I times 
by parts then using Laplace’s method. 

[82] The product of the four unitary operators under the trace 
represents in the single-particle Hilbert space a rotation 
of angle <5 around some axis. It is easy to explicitly eval¬ 


uate this trace as a function of the angles 9,9',4>,4>' in 
the case 1=1, where each operator can be replaced by 
a well-known 3x3 rotation matrix in the usual, three- 
dimensional space. This leads to the expression (IB57II . 

[83] Furthermore the Hilbert space was limited in reference 
[l^ to the kets \x,u, i,mz) with a: > 0, as the sym¬ 
metry condition (IB58II allows, which amounts to adding 
an extra projector Px>o- This complicates things be¬ 
cause Px>o and U do not commute. Fortunately, in 
calculating operator traces, one can use the properties 
UP = PU = P, VxKO = U'Px>oU and Vxko + Px>o = 1, 
as well as the invariance of the trace under a cyclic 
permutation of the operators, so that Tr('Pa,<oPAP) = 
Tr{UPx>oUPAP) = Tr{Vx>oPAP) = i Tr(PAP) = 
iTr(AP) and Tr(APP,<oPP) = Tr(APP,>oPP) = 
i Tr(APP), where A and B are arbitrary operators. 

[84] If one writes the operator M of reference [1^ before its 
restriction to the subspace of symmetry (IB58I) as P -I- 
Ko + UKoU then after implementation of the symmetry 
and restriction to the Hilbert space of kets \x,u,£,mz) 
with a; > 0, it becomes Pxyol'P + {1 + U)Ko{l + U)]Px>o- 
The first, second and third contributions in the right- 
hand side of Eq. (IB60I) are respectively given by —{£ + 
1/2)/{2tv) times the integral over S G R of Tti^e{KoU), 
of -^Tie,,[Ko{UKoU)] and ol -Tte,e{KoU), where we 
have set Kq = T>~^Ko and the index £,£ means that the 
trace is restricted to the states \£,m) with (—1)"“ = e. 
Note that UPU = T) and [Pa:>o,P] = 0. 

[85] The £ independent function in factor of the sine function 
in the integrand of Eq. (IB62I) is a smooth function of 5 
over ]0, tt] , with all its even order derivatives (including 
the zeroth order) vanishing at <5 = 0 and all its odd order 
derivatives vanishing at <5 = tt. Under repeated integra¬ 
tion by parts (always integrating the sine function), the 
fully-integrated term vanishes at the boundaries and one 
pulls out at each step a factor (£ -|- 1/2)“^. So Eq. ([B62I) 
is 0[{£ + 1/2)“"] when £ —^ -|-oo, for all integers n. 


















